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PREFACE 

This  volume  is  a  summary  of  presentations  and  dis¬ 
cussions  of  a  technical  workshop  on  Radar  Scattering  from 
Random  Media,  held  at  the  Institute  for  Pure  and  Applied 
Sciences,  University  of  California  (San  Diego),  La  Jolla, 
California,  on  5  -  16  August  1968,  and  sponsored  by  the 
Advanced  Research  Projects  Agency.  The  Workshop  was 
divided  into  Theoretical  and  Experimental  Panels .  Summaries 
of  the  reports  of  these  Panels  are  the  result  of  collaboration 
among  several  Workshop  participants.  Special  thanks  are  given 
to  Dr.  S.  Rand  (IPAS)  and  to  Dr.  S.  C.  Lin  (IPAS)  for  preparing 
and  editing  the  panel  reports.  A  complete  transcript  of  the 
proceedings  of  the  Workshop  was  prepared  by  Dr.  R.  Ruquist 
(MIT /LL) ,  and,  although  not  included  herein,  is  available 


upon  request.  Special  thanks  is  given  to  Dr.  T.  0.  Philips 
(BTL),  general  editor  of  this  report,  who  skillfully  molded 


the  diverse  sections  of  the  report  into  a  coherent  document 


K.  Kresa  (ARPA) 


Workshop  Chairman 


ii 


CONTENTS 

Page 


Preface  i 

Table  of  Contents  11 

Chapter  1.  Introduction  1 

A.  The  Scope  of  the  Workshop  and  of  this  report  1 

B.  Major  Conclusions  5 

Chapter  2.  Workshop  Participants  and  Program  8 

Chapter  3*  Theoretical  Studies  of  Scattering  from  12 

Random  Media 

A.  Introduction  12 

B.  Perturbation  Techniques  18 

C.  Watson's  Transport  Equation  31 

D.  Information  Theory  Formulation  35 

E.  Computer  Experiments  39 

Chapter  4.  Laboratory  Facilities  for  Scattering  44 

Experiments 

f 

Chapter  5-  Conclusions  and  Recommendations  49 

A.  Conclusions  and  Recommendations  of  the  49 

Theoretical  Panel 

B.  Conclusions  and  Recommendations  of  the  54 

Experimental  Panel 


Appendix  A.  Bibliography  on  Theory  and  Experiments  in 
Scattering  from  Turbulent  Plasmas. 
Appendix  B.  Small-scale  Structure  and  Viscous  Cutoff 
in  Scalar  Spectrum  of  Hypersonic  Wake 
Turbulence  (by  Shao-Chi  Lin) 

Appendix  C.  The  Compatibility  of  Electromagnetic 
Scattering  Theory  and  Field  Data 
(SECRET;  bound  and  distributed 
separately)  . 


I 


- 1  - 


Chapter  1 
INTRODUCTION 

| 

Section  1.  A  The  Scope  of  the  Workshop  and  of  This  Report 
The  understanding  of  radar  returns  from  turbulent 
plasmas  has  been  the  focus  of  much  interest  for  many  years.  This 
subject  is  important  because  radar  is  extensively  utilized  in  both  the 
field  and  the  laboratory  as  a  remote  sensor  to  provide  basic 
information  on  turbulent  flows  and  many  aspects  of  the  physics  of 
reentry. 

The  use  of  radar  scattering  in  this  diagnostic  capacity  is 
difficult,  since  radar  scattering  provides  only  an  indirect  method 
of  probing  the  fluid  dynamics  and  turbulent  structure  of  the 
medium.  *a  a  weakly-ionized  turbulent  plasma  (such  as  a  turbulent 
wake)  microwave  radiation  is  scattered  by  the  free  electrons  and 
one  must  know  or  guess  the  relationship  between  the  electronic 
phenomena  observed  and  the  underlying  fluid  dynamics  and  turbulent 
structure.  Moreover,  in  order  even  to  interpret  radar  scattering 
measurements  in  terms  of  electron  densities  and  electron  motions, 
it  is  usually  necessary  to  have  some  theoretical  understanding  of  the 


electromagnetic  scattering  process. 
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In  certain  limiting  cases  (for  example,  in  the  case  of 
low  mean  electron  densities  and  small  electron  density 
fluctuations)  simple  and  reliable  theories  are  available.  But  it 
is  seldom  easy  to  determine  the  conditions  under  which  these 
simple  theories  should  be  valid  or  the  nature  of  the  deviations  from 
those  theories.  Other  theories  which  may  be  more  widely  applicable 
are  frequently  so  complex  that  they  are  not  useful  when  detailed 
quantitative  results  are  required. 

This  Workshop  was  motivated  by  ARPA's  desire  to 
assess  the  present  state  of  knowledge  of  scattering  theory  as 
applied  primarily  to  radar  scattering  from  turbulent  wakes. 

Both  theoretical  and  laboratory  experimental  research  were 
reviewed  and  pertinent  correlations  between  theory  and 
experiment  were  examined.  In  addition,  the  results  obtained 
from  radar  measurements  were  presented  and  the  relevance  of 
present  theories  to  the  understanding  of  these  data  were  discussed. 

The  goal  of  the  Workshop,  and  the  purpose  of  this 
report,  was  to  prepare  a  critical  review  of  the  present  research 
in  this  area  with  respect  to  the  understanding  of  radar 
scattering  from  ionized  turbulent  media  and,  in  addition,  to 
recommend  areas  for  future  emphasis  and  support. 
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The  activities  of  the  Workshop  were  highly  informal, 
as  is  this  volume  of  reports  emanating  from  the  Workshop. 
Although  several  participants  were  invited  to  prepare  lectures 
on  relevant  topics,  most  of  the  Workshop  sessions  consisted 
of  informal  discussion  and  debate.  There  was  no  attempt  to 
record  these  discussions  in  detail  herein,  and  this  present 
volume  i s  not  intended  to  be  a  thorough  documentation  of  the 
material  considered  by  the  Workshop.  Instead  it  should  be 
viewed  as  a  collection  of  written  reports,  prepared  by 
individuals  or  groups  present  at  the  Workshop,  which  may  be 
of  some  interest  or  use  to  a  wider  audience. 

Section  B  of  this  introductory  Chapter  1  summarizes 
some  particularly  significant  conclusions  and  recommendations 
of  the  Workshop.  A  list  of  the  Workshop  participants  and  a 
summary  of  the  program  are  found  in  Chapter  2. 

Chapter  3  reviews  the  present  status  of  the  theory 
of  incoherent  electromagnetic  scattering  from  random  media. 

The  theoretical  models  reviewed  are  those  which  were  dis¬ 
cussed  at  the  Workshop;  those  which  appear  most  useful  for 
the  understanding  of  wake  scattering  received  the  most  attention 
both  in  the  Workshop  discussions  and  in  this  report.  Table 
I  (in  Chapter  3)  presents  a  useful  summary  and  comparison  of 
these  various  theories.  Chapter  4  contains  a  brief  description 
of  some  of  the  laboratory  facilities  in  which  scattering  experi¬ 
ments  are  being  performed.  Results  from  these  scattering 
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experiments  were  presented  to  the  Workshop  but,  since  they 
are  well  documented  in  the  technical  literature,  these 
results  are  not  included  here. 

During  the  last  three  days  the  Workshop  was 
divided  into  a  Theoretical  Panel  and  an  Experimental  Panel, 
according  to  the  research  interests  of  the  participants. 

Each  panel  reviewed  the  Workshop  discussions  and  prepared  a 
set  of  conclusions  and  recommendations  which  are  to  be 
found  In  Chapter  5.. 

Three  appendices  are  attached  to  this  report. 
Appendix  A  contains  a  bibliography  of  papers  on  the  theory 
of  scattering  from  random  media,  as  well  as  papers  document¬ 
ing  the  results  of  laboratory  scattering  experiments.  It  is 
based  upon  a  bibliography  of  theoretical  papers  prepared  for 
the  Workshop  by  A.  Hochstim;  additional  entries  were  supplied 
by  T.  0.  Philips.  Appendix  B  is  a  paper  by  S.  C.  Lin  on  the 
spectrum  of  electron  density  fluctuations  in  hypersonic  wakes. 
It  is  an  expanded  version  of  a  lecture  given  at  ;he  Workshop. 
Appendix  C  represents  an  attempt  to  determine  the  applica¬ 
bility  of  various  theories  of  scattering  from  random  media  to 
the  interpretation  of  field  data.  This  paper  is  the  output  of 
a  "working  session"  of  the  Workshop,  which  was  organized  by 
S.  Edelberg,  Some  of  the  conclusions  reported  in  Section  B  of 
Chapter  5  are  based  upon  this  3tudy.  Appendix  C  is  the  only 
classified  section  of  this  Workshop  Report  and  is  therefore 
being  distributed  separately. 
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Section  1.  B  Major  Conclusions 

Chapter  5  contains  conclusions  and  recommendations 
which  were  identified  by  the  Theoretical  and  Experimental 
Panels  of  the  Workshop.  A  few  of  these  conclusions  appear  to 
be  especially  significant  for  future  investigations  of  scattering 
from  reentry  wakes. 

Validity  of  the  first  Born  approximation.  The  simplest 
and  most  widely  used  theory  of  electromagnetic  scattering  from 
random  media  is  the  first  Born  approximation;  the  conditions  for 
the  validity  of  this  theory  have  been  investigated  by  Salpcter  and 
Treiman.  Several  different  sources  reviewed  at  the  Workshop 
suggest  that  these  Salpeter-Trieman  conditions  may  be  unnecessarily 
restrictive  for  the  purposes  of  scattering  studies  where  agreement 
of  theory  and  experiment  within  a  factor  of  two  is  adequate. 

For  example,  the  computer  experiments  of  Hochstim  and  the  laboratory 
experiments  of  Guthart  and  his  colleagues  suggest  that  the  first 
Born  approximation  may  be  valid  (to  within  a  factor  of  two)  for 
electron  density  fluctuation  levels  which  approach  the  critical 
electron  density.  (On  the  other  hand,  the  experiments  of 
Granatstein  suggest  that  significant  deviations  from  the  first  Born 
approximation  may  occur  at  somewhat  lower  levels  of  electron  density 


fluctuations. ) 
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A  brief  study  performed  at  the  Workshop,  using  radar  wake 
scat*-  lng  data  from  reentry  experiments,  and  using  a 
theoretical  approach  due  to  Shkarofsky  and  to  Feinstein, 
also  indicates  that  the  first  Born  approximation  may  be  valid 
under  a  wider  range  of  reentry  conditions  than  had  previously 
been  realized.  It  will  require  further  theoretical  and 
experimental  studies  to  verify  that  the  first  Born  approxi¬ 
mation  is  indeed  so  widely  applicable  to  the  wake  scattering 
problem. 

Use  of  the  distorted  wave  Born  approximations. 

There  are  several  generalizations  of  the  first  Born  approxi¬ 
mation  which  have  appeared  in  the  scattering  theory  literature. 
Most  of  these  generalizations  are  equivalent  to  one  of  the 
five  "distorted  wave  Born  approximations"  which  are  discussed  ' 
in  Chapter  3>  where  they  are  given  the  names  DWBA-1,  -2a,  -2b, 
-2c,  and  the  "Kraichnan  model".  Each  of  these  theories  is  a 
special  case  of  the  theories  which  follow  it  in  the  list.  It 
was  concluded  by  the  Workshop  participants  that  there  may  be 
little  point  in  applying  the  theories  DWBA-1,  -2a,  or  -2b  to 
the  wake  scattering  problem,  since  it  appears  that  the  model 
DWBA-2c  is  only  slightly  more  complicated  and  would  be 
expected  to  yield  significantly  better  quantitative  results. 
Other  theories,  such  as  the  Kraichnan  model  and  the  Information 
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Theory  Formulation  (see  Section  2.D)  may  in  principle  be 
more  general,  but  do  not  appear  to  yield  practical  results 
for  the  study  of  wake  scattering. 

Definition  of  the  scattering  medium.  Even  if  one 
has  a  collection  of  theoretical  models  for  scattering  from 
a  random  medium,  accompanied  by  conditions  for  the  validity 
of  these  models,  the  choice  of  the  most  useful  model  depends 
upon  the  nature  of  the  scattering  medium.  Is  the  medium, 
for  example,  adequately  characterized  by  a  statistical  des¬ 
cription  in  terms  of  correlation  functions  of  random  variables? 

Is  it  necessary  to  specify  in  detail  the  nature  of  the  boun¬ 
dary  of  the  medium?  Are  the  gradients  of  mean  electron 
density  sufficiently  small  that  diffraction  effects  may  be 
neglected? 

It  is  not  clear  that  the  currently  used  statistical 
descriptions  of  a  turbulent  medium  are  always  adequate;  under 
some  conditions  it  may  be  necessary  to  have  a  more  complete 
description  of  a  single  realization  of  a  random  medium.  The 
development  of  alternative  characterizations  of  the  scattering 
medium  should  be  a  major  goal  of  future  studies,  both  theoretical 
and  experimental,  of  the  fluid  dynamics,  the  turbulence  proper¬ 
ties,  and  the  chemistry  of  wakes. 
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Chapter  2 

WORKSHOP  PARTICIPANTS  AND  PROGRAM 


The  following  persons  participated  in  one  or  more 

of  the  sessions  of  the  Workshop: 

KEITH  A.  BRUECKNER,  University  of  California,  San  Diego 
ANTHONY  DEMETRIADES,  Aeronutronic  Division,  Philco 

Corporation 

SEYMOUR  EDELBERG,  Lincoln  Laboratory,  Massachusetts 

Institute  of  Technology 

LEOPOLD  B.  FELSEN,  Polytechnic  Institute  of  Brooklyn 
CARL  H.  GIBSON,  University  of  California,  San  Diego 
VICTOR  L.  GRANATSTEIN,  Bell  Telephone  Laboratories 
HAROLD  GUTHART,  Stanford  Research  Institute 
ADOLF  HOCHSTIM,  Institute  for  Defense  Analyses 
JOHN  JAREM,  Drexel  Institute  of  Technology 
TUDOR  W.  JOHNSTON,  RCA  Victor  Research  Laboratory 

(Montreal) 

KENT  KRESA,  Advanced  Research  Projects  Agency 
S.  C.  LIN,  University  of  California,  San  Diego 
L.  R.  MARTIN,  Lincoln  Laboratory,  Massachusetts  Institute 
of  Technology 

RAYMOND  F.  MISSERT,  Cornell  Aeronautical  Laboratory 
MARVIN  H.  MITTLEMAN,  University  of  California,  Berkeley 
S.  S.  PENNER,  University  of  California,  San  Diego 
THOMAS  0.  PHILIPS,  Bell  Telephone  Laboratories 
ROBIN  I.  PRIMICH,  A.C.  Electronics  Defense  Research 

Laboratories 

ANDREW  PROUDIAN,  Heliodyne  Corporation 
S.  RAND,  University  of  California,  San  Diego 
JACQUES  RENAU,  Aerospace  Corporation 

RICHARD  S.  RUFFINE,  Advanced  Ballistic  Missile  Defense 

Agency 

RICHARD  RUQUIST ,  Lincoln  Laboratory,  Massachusetts 

Institute  of  Technology 

I.  P.  SHKAROFSKY,  RCA  Victor  Research  Laboratory 

(Montreal) 

BURTON  STROM,  Riverside  Research  Institute 
K.  SULZMANN,  University  of  California,  San  Diego 
EMMETT  A.  SUTTON,  Avco-Everett  Research  Laboratory 
KENNETH  M.  WATSON,  University  of  California,  San  Diego 
S.  ZIVANOVIC,  A.C.  Electronics  Defense  Research 
Laboratories 
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The  Workshop  sessions  included  formal  lectures, 
spontaneous  presentations,  informal  discussions,  and  working 
sessions  (which  sometimes  involved  the  entire  Workshop  and 
sometimes  only  a  subcommittee  of  the  attendees).  In  the 
following  outline  we  attempt  to  refer  to  all  of  the  signifi¬ 
cant  discussions  (and  not  only  to  the  more  formal  lectures). 
We  have  assigned  titles  to  the  informal  and  spontaneous  dis¬ 
cussions;  the  titles  will,  we  hope,  characterize  the  nature 
of  the  subjects  discussed. 

Monday,  August  5,  1968 
Registration. 

Goals  for  the  Workshop,  K.  Kresa. 

Definitions,  posing  of  the  problem,  and  review  of  scattering 
theory  literature.  A,  Hochstim. 

Review  of  scattering  from  random  media,  J.  Keller  (paper 
presented  by  J.  Jarem) . 

Discussions : 

Comments  on  a  paper  by  Frisch,  J.  Jarem. 

The  Tatarskii  approach  to  scattering  theory,  L.  Felsen. 
Identification  of  physical  parameters  pertinent  for 
scattering  phenomena,  R.  Ruffine  and  others. 
Comments  on  the  conditions  for  various  scattering 
theories,  I.  Shkarofsky. 

Tuesday,  August  6 

Derivation  of  the  transport  equation  for  scattering  from  a 
random  medium,  K.  Watson. 

A  one-dimensional  scattering  experiment  on  the  computer, 

A.  Hochstim. 

Wednesday ,  August  7 

A  hueristic  multiple  scattering  theory,  I.  Shkarofsky. 

The  distorted  wave  Born  approximation,  M.  Mittleman. 

A  solution  for  scattering  from  a  one-dimensional  slab, 

J.  Jarem. 
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Working  session:  Numerical  values  of  physical  parameters 
pertinent  for  scattering  phenomena  (R.  Rufflne  and 
others ) . 

Discussions : 

The  applicability  of  scattering  theories  to 
scattering  from  reentry  vehicle  wakes. 

Summary  of  available  theoretical  methods,  A.  Hochstlm. 

A  theorem  of  Lax,  K.  Watson. 

Distorted  wave  Born  approximations,  J.  Jarem,  L.  Felsen. 
Limitations  to  transport  theory,  K.  Watson. 

Thursday,  August  8 

Distorted  wave  Born  approximations,  L.  Felsen. 

Statistical  methods  In  scattering  theory,  A.  Proudlan. 
Properties  of  wakes  of  hypersonic  projectiles,  E.  Sutton. 

(with  comments  by  A.  Demetrlades ) . 

Discussions:  Necessary  extensions  to  current  theory. 

Working  session;  Preparation  of  reports  on  the  status  of 
theories  of  scattering  from  turbulent  plasmas. 

Friday,  August  9 

Working  session:  Preparation  of  reports  on  the  status  of 
theories  of  scattering  from  turbulent  plasmas. 

Monday,  August  12 

Summary  of  theories  of  scattering  from  turbulent  plasmas, 

A.  Proudlan. 

Discussions:  The  current  status  of  scattering  theory. 

Experiments  on  scattering  from  a  turbulent  plasma  jet, 

H.  Guthart . 

Characteristics  of  turbulent  wakes,  A.  Demetrlades. 

Diagnostic  study  of  a  turbulent  plasma  jet,  T.  Johnston. 
Scattering  experiment  with  a  flowing  turbulent  plasma, 

V.  Granatsteln. 

Tuesday,  August  13 

The  spectrum  of  wake  turbulence  -  Comparisons  of  field  data 
with  theory,  E.  Sutton. 

Recent  TRADEX  wake  scattering  observations,  L.  Martin. 

Recent  wake  scattering  observations  at  L-band  and  C-band, 

R.  Missert. 

Discussions : 

A  multiple  scattering  criterion  from  Ionospheric  physics, 
J.  Renau. 

The  limits  of  the  Born  approximation  based  upon  cross- 
polarized  scattering,  R.  Rufflne. 
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Wednesday,  August  Ik 

Working  session  (Experimental  panel):  Comparison  ol  theoret 
leal  scattering  models  with  recent  field  data, 

S.  Edelberg  and  others. 

Thursday,  August  15 

The  spectrum  of  electron  density  fluctuations  In  turbulent 
reentry  wakes,  S.  Lin. 

Working  session  (Theoretical  and  experimental  panels): 
Formulation  of  the  conclusions  of  the  Workshop  and 
recommendations  for  future  research. 
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Chapter  3 

THEORETICAL  STUDIES  OF  SCATTERING 
FROM  RANDOM  MEDIA 

Section^  .A  Introduction 

During  the  first  four  days  (August  5  through  8)  the 
Workshop  heard  lectures  on  a  number  of  different  approaches 
to  the  theory  of  scattering  from  random  media.  The  principal 
goals  during  these  four  days  were  (1)  to  achieve  an  under¬ 
standing  of  the  relationships  among  the  various  approaches, 

(II)  to  estimate  the  range  of  validity  of  each  approach,  and 

(III)  to  assess  the  likelihood  that  any  given  approach  could 
be  usefully  applied  to  the  specific  problem  of  wake  scattering. 
Much  of  the  information  assembled  concerning  these  questions 

is  presented  in  Table  I.  A  summary  of  the  general  conclusions 
and  recommendations  concerning  theoretical  studies  is  given  in 
Chapter  5. 

In  most,  but  not  all,  cases,  the  various ‘approaches 
to  scattering  calculations  were  discussed  in  the  context  of  a 
scalar  wave  equation.  Little  detailed  study  was  given  to  the 
possibility  of  extending  the  techniques  to  electromagnetic 
scattering  problems  in  which  polarization  effects  are  impor¬ 
tant;  however,  some  assessment  of  the  possibility  of  such  an 
extension  was  included  in  Table  I. 
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TABLE  I.  SCATTERING  THEORY  STATUS 


Range  of  Validity 

Application  to 
Slab 

Finite 

Geometry 

First  Born 

«  i 
e 

ObD  «  1 

yes 

yes 

DWBA-1 

—  «  1 
e 

aBL  «  l  afiD  «  l 

done  1-D  Jarem  3D 
in  progress 

* 

hard  for  exact 
Green's  function 

Heuristic  Model 
(DWBA-2a) 

Ss.  «  1 

e 

yes 

yes  -  Shkarofsky 

Simplified  Model 
(DWBA-2b) 

R  «  1  f 

half-plane  in 
process  -  Felsen 

probably  with 
much  work 

Full  Model 
( DWBA-2c ) 

R  «  1  f 

1-D  possibly 
Mittleman 

* 

very  hard 

Transport  Theory 

kD  »  1 

R  «  1 

V  In  n  «  k 

has  been  done  in 
radiative  transfer 

— > 

Kraichnan 

R  unbounded 

* 

can  be  done  -  hard 

* 

hard 

Information  Theory 
Formulation 

R  unbounded 

Definitions : 

R  ■=  k  ij  L  *  correlation  length  k  =  2w/X 

D  =  physical  dimension  along  beam 

=  total  Born  cross  section  per  unit  volume 

n 


Ray  tracing  (V  In  n  «  k)  makes  problem  more  tractable. 

^  Watson  claims  models  give  correct  backscatter  power  to  a  factor  of  2 
provided  R  «  1  for  simplified  model  and  y/R  <  1  for  full  model. 
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TABLE  I.  SCATTERING  THEORY  STATUS  (Cont'd) 


"  ' 

Status  of 
Research 

Verification 

Polarization 

Effects 

First  Bor;i 

complete 

with  experiment 

yes  -  none  for 
k 

first  Born;  6 
effect  2nd  Born 

DWBA-1 

slab  with  3-D 
formulated  (Jarem) 

1-D  Hochstim 

3-D  none 

hasn't  been  done, 
but  possible 

Heuristic  Jtodel 
(DWBA-2a) 

extend  to  include 
density  profiles, 
anisotropy 

pending 

yes 

Simplified  Model 
(DWBA-2b) 

? 

? 

from  curved  paths 
only  or  2nd  Bom 

Full  Model 
(DWBA-2c) 

formulated 

1-D  underway 
Felsen,  Mittleman 

from  curved  paths 
only  or  2nd  Bom 

Transport  Theory 

diffusion  limit 
done  -  ready  for 
add'l  applications 

probably  in  other 
problems  -  diffusion 
limit  will  be  check¬ 
ed  during  workshop 

yes 

Kraichnan 

formulated 

verified  on 
other  models 

yes 

Information  Theory 
Formulacion 

formulated 

? 

should  give 
correct  value 
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It  was  decided  to  defer  to  some  other  occasion  any 
discussion  of  the  use  of  the  various  scattering  theories  to 
calculate  the  frequency  spectrum  of  waves  scattered  from  ran¬ 
dom  media.  A  discussion  of  frequency  effects  is  important 
and  should  be  attempted  as  a  sequel  to  this  Workshop. 

In  order  to  attempt  to  understand  the  ensemble- 
averaged  microwave  scattering  from  a  turbulent  wake,  four 
general  approaches  were  considered:  (i)  perturbation  tech¬ 
niques,  (ii)  Watson's  transport  equation,  (iii)  the  so-called 
information  theory  formulation,  and  (iv)  computer  experiments. 

The  approach  which  received  the  most  attention  in¬ 
volved  a  variety  of  modifications  of  the  perturbation  scheme, 
whereby  an  expansion  is  performed  in  terms  of  a  small  param¬ 
eter.  The  prototype  method,  beyond  which  generalizations 
were  proposed,  is  the  well-known  "first  Born  approximation," 
an  attractive  starting  point  because  of  its  simplicity  as  well 
as  its  applicability  to  a  sufficiently  underdense  wake. 
Attempts  to  generalize  the  first  Born  approximation  included 
a  number  of  techniques  which  are  collectively  referred  to  as 
"distorted  wave  Born  approximations";  in  these  techniques  the 
expansion  parameter  is  the  electron  density  fluctuation, 
rather  than  tne  total  electron  density  as  in  the  first  Born 
approximation.  These  methods  are  a  useful  improvement  when 
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there  are  regions  of  the  wake  where  the  mean  electron  density 
approaches,  or  possibly  even  exceeds,  the  critical  value 
corresponding  to  the  Incident  microwave  frequency,  while  the 
statistical  fluctuations  from  the  mean  are  much  less  than  the 
critical  density.  These  perturbation  methods  will  be  dis¬ 
cussed  In  Section  3.B. 

A  second  approach  to  the  scattering  problem  was  out¬ 
lined  In  a  lecture  by  K.  M.  Watson  In  which  he  discussed  the 
derivation  of  a  vector  transport  equation  for  the  Incoherently 
scattered  energy.  The  derivation  starts  from  the  multiple 
scattering  equations  for  scattering  from  a  collection  of  N 
Individual  electrons.  This  approach  has  the  merit  that  It 
could  take  advantage  of  the  many  existing  solutions  to  the 
transport  equation.  Watson’s  equation  is  discussed  In 
Section  3.C. 

Limited  attention  was  given  to  a  completely  differ¬ 
ent  approach,  which  has  been  called  the  "Information  Theory 
Formulation,"  a  technique  which  Is  of  the  nature  of  a  vari¬ 
ational  method.  The  theory  as  presented  to  the  Workshop  Is 
in  a  primitive  stage,  but  Is  outlined  In  Section  3. D. because 
it  has  a  slim  possibility  of  being  the  only  true  theory  which 
may  eventually  be  used  to  describe  scattering  from  a  wake 
with  overdense  fluctuations.  The  discussion  was  brief  because, 
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when  applied  to  the  wave  equation  with  a  random  variable, 
the  formulation  has  not  reached  the  point  where  much  can  be 
surmised  about  its  applicability. 

A  fourth  class  of  techniques  was  proposed  which 
included  either  numerical  methods,  that  is,  "computer  experi¬ 
ments,"  or  actual  scattering  experiments  to  be  performed  on 
models.  None  of  the  models  would  simulate  the  wake,  but 
rather  they  would  be  simple  mocked-up  situations,  with  sta¬ 
tistics  which  are  completely  prescribed,  in  order  to  test  the 
various  scattering  theories.  The  advantage  is  clear  in  that 
the  complicated  wake  structure  is  ignored,  so  that  the  testing 
is  limited  entirely  to  the  scattering  theories.  The  various 
methods  for  extending  some  simple  one-dimensional  computer 
experiments,  which  have  already  been  performed,  will  be 
described  in  Section  3*E. 

An  obvious  omission  in  the  discussions  at  this 
Workshop,  by  consensus  of  the  participants,  involved  theories 
of  scattering  from  overdense  random  surfaces.  It  is  hoped 
that  these  theories,  which  are  expected  to  be  appropriate  to 
wake  scattering  under  some  reentry  conditions,  will  be  the 
subject  of  another  meeting. 
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Section  3.B  Perturbation  Techniques 

Five  proposals  were  considered  for  extending  the 

# 

first  Born  approximation.  The  first  Born  approximation 
assumes  that  the  coherent  wave  propagates  In  the  random 
medium  with  the  free  space  wave  number  k  .  (The  coherent 
wave  is  then  exactly  the  wave  Incident  upon  the  medium.) 

In  all  of  the  proposed  extensions  to  the  theory  the  coherent 
wave  propagates  with  an  effective  wave  number  k(x)  which 
differs  from  the  free  space  value  kQ.  These  extensions  are 
all  referred  to  as  "distorted  wave  Born  approximations." 

In  the  first  such  model,  which  will  be  referred  to 
as  the  "simple  model"  (or,  in  the  language  used  by  the  Work¬ 
shop,  DWBA-1),  the  effective  wave  number  k(x)  differs  from 
kQ  because  of  the  mean  electron  density.  The  effects  of 
electron  density  fluctuations  on  the  coherent  field  are 
ignored . 


Some  calculations  have  been  performed  in  the  second  Born 
approximation  (that  is,  second  order  terms  of  the  Born 
perturbation  series  are  retained).  However,  even  these 
calculations  are  considerably  more  difficult  than  the  first 
Born  approximation  and,  except  in  very  special  cases,  it 
seems  hopeless  to  expect  to  be  able  to  calculate  hi  'her 
order  terms  in  the  Born  series. 
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The  other  four  "distorted  wave  Born  approximations" 
include  some  of  the  effects  of  electron  density  fluctuations 
on  the  coherent  wave.  The  "heuristic  model"  (DWBA-2a)  in¬ 
cludes  in  an  ad  hoc  manner  an  attenuation  of  the  v.  oherent 
wave;  this  attenuation  is  the  result  of  the  scattering  of 
energy  out  of  the  coherent  beam  by  electron  density  fluctu¬ 
ations.  The  fluctuations  therefore  modify  the  imaginary  part 
of  the  effective  wave  number  k(x).  In  the  "expanded  model" 
(DWBA-2b)  the  real  part,  as  well  as  the  imaginary  part,  of 
the  effective  wave  number  is  modified  by  fluctuations  to 
first  order  in  the  scattering  coefficient.  The  "full  model" 
(DWBA-2c )  brings  in  the  effects  of  higher  order  scatters  in 
the  equation  for  the  coherent  field,  but  ignores  them  in  the 
equation  by  which  the  Incoherent  field  is  derived  from  the 
coherent  field.  Although  there  appears  to  be  an  inconsist¬ 
ency  in  this  model,  in  that  terms  are  ignored  which  seem  to 
be  as  large  as  terms  which  are  retained,  this  model  is  taken 
seriously  not  only  because  of  the  reasonableness  of  the 
formulation,  but  also  because  it  has  been  shown  to  result  in 
improvement  over  the  simplified  theories  when  applied  to 
nuclear  scattering  problems.  Finally,  the  "Kraichnan  model" 
includes  additional  multiple  scatters  beyond  the  full  model, 
and  appears  to  be  somewhat  of  an  improvement.  However,  it 


20 


involves  a  nonlinear  integral  equation,  which  in  most  appli¬ 
cations  would  be  extremely  difficult  to  solve. 

The  five  perturbation  schemes  described  above  in¬ 
volve  increasing  complications,  and  each  includes  the  pre¬ 
vious  models  as  special  cases.  A  rather  direct  derivation 
of  the  "full  model"  (DWBA-2c)  is  available,  and  this  will  be 
outlined  here. 

Consider  the  equation 

(Lq+L1  )ij>  =  0  ,  (1) 

where  Lq  and  are  linear  operators,  Lq  being  sure  and 
being  stochastic.  The  electron  density  fluctuations  are 
considered  to  be  a  random  process  and  are  included  in  L^. 

The  mean  electron  density  may  appear  either  in  (as  in  the 
Born  approximation)  or  in  Lq  (as  in  the  various  distorted 
wave  Born  approximations).  Because  of  the  total  field 
will  contain  a  fluctuating  (stochastic)  part  5^  as  well  as 
a  coherent  part  4> , 


ip  =  <p  +  Sip  .  (2) 

By  substituting  equation  (2)  into  equation  (1), 


and  averaging  over  the  random  variables,  one  obtains 
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Lo4>  +  <L  =0  (3) 

where  it  will  be  assumed  that  is  centered,  that  Is, 

<(L^>  =  0,  and  therefore  also  =  0.  By  subtracting 

equation  (3)  from  equation  (1),  we  have  the  equation  for 
associated  with  a  given  realization, 

Lq6^  +  L14»  +  -  <L1«4»>  =0.  (H) 

The  fundamental  assumption  of  DWBA-2c  is  that  the  quantity 
1^64*  -  Is  much  smaller  than  the  other  terms  of  equa¬ 

tion  (4),  and  may  therefore  be  discarded.  This  Is  an  ad  hoc 
assumption,  with  no  theoretical  justification  offered  by  any 
of  the  participants  at  the  Workshop.  It  is  accepted  there¬ 
fore,  that  the  "full  model"  DWBA-2c  may  constitute  no  Im¬ 
provement  over  the  "simple  model"  DWBA-1  In  which  additional 
terms,  which  appear  to  be  of  the  same  order  of  magnitude,  are 
neglected.  Although  the  a  priori  expansion  parameter  of  both 
DWBA-1  and  DWBA-2c  (as  well  as  of  the  "expanded  model" 
DWBA-2b,  of  intermediate  complexity)  would  appear  to  be  the 
same,  namely  the  electron  density  fluctuations,  the  model 
DWBA-2c  Is  expected  to  be  an  Improvement  because  it  has  given 
results  which  compare  better  with  experiment  In  the  corre¬ 
sponding  nuclear  scattering  problem.  Thus  for  DWBA-2c,  equa¬ 
tion  (4)  is  replaced  by 
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Lo5if<  +  L14>  =  0  (5) 

with  the  formal  solution, 

S\p  =  -L"1  L^.  (6) 

When  equation  (6)  is  combined  with  equation  (3),  we  have  a 
single  equation  for  the  coherent  field, 

[lq  -  <LXL^  Lx>]  4>  =0.  (7) 

The  quantity  <L1L~1  L^)>  is  an  integral  operator  which  in¬ 
volves  the  statistical  properties  of  the  electron  density 
fluctuations.  This  term  is  ignored  in  both  the  first  Born 
approximation  and  the  "simple"  distorted  wave  Born  approxi¬ 
mation  (DWBA-1). 

In  the  first  Born  approximation  the  sure  operator 
L  is  taken  to  be 


L  = 

o 


(8a) 


where  kQ  is  the  free  space  wave  number  (kQ 


2TTfQ/c)  of  the 


monochromatic  plane  wave  incident  upon  the  random  medium; 
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the  wave  equation  for  the  coherent  field  $  is 


J  =  {v2+kl)*  =  °* 


(8b) 


In  the  model  DWBA-1,  the  incident  free  space  wave 
number  kQ  is  replaced  by  an  effective  wave  number  k(x),  where 
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(9a) 


which  includes  the  effects  on  the  coherent  wave  of  the  mean 
electron  density  <n(x)>.  Here  n  is  the  critical  electron 
density  corresponding  to  the  incident  frequency  fQ.  Then 
the  wave  equation  for  the  coherent  field  is 


V  -  [v2  +  k(x  ) 2]  4)  =  0.  (9b) 

It  is  known  that  equation  (9b)  can  be  solved  for  the  coherent 

field  in  a  number  of  simple  geometries,  even  when  4>  is  a 

vector  field.  The  problem  of  an  infinite  cylinder,  with 

k(x)  a  step  function,  is  exactly  solvable. 

A  major  problem  of  applying  the  full  model,  given 

by  equation  (7),  is  associated  with  the  determination  of  L-1, 

o 

that  is,  of  obtaining  the  Green's  function  for  a  real  problem. 
For  example,  if  one  writes 
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Lx  =  k(x)2p(x) 


(10) 


where  p(x)  Is  a  stochastic  variable  which  is  proportional 
to  the  electron  density  fluctuation,  then  equation  (7)  may 
be  written 


0  *  [v2  +  k(x)2]$(x; 


k(x)2k(x'  )2G^°^  (x,x'  )<(y(x)y(x'  )}>4>(xf  )d^xf  (11) 


where  L-1  =  G^0^(x,x’)  is  the  Green’s  function.  Even  if  k( 
o 

is  constant  throughout  the  wake,  the  Green’s  function  is 
extremely  complicated  (although  available  for  the  case  of  a 
cylinder)  because  of  the  finite  geometry.  It  was  proposed 
at  the  Workshop  therefore,  to  use  the  Green’s  function 
associated  with  an  infinite  space,  but  with  wave  number  k 
which  corresponds  to  some  mean  electron  density  inside  the 
wake.  The  assumption  is  that  the  correlation  length  and  the 
dimension  of  the  wake  are  both  large  compared  to  the  radia¬ 
tion  wavelength  inside  the  wake.  Using  the  proposed  assump¬ 
tion,  we  have 


G(o)(x,x’) 


exp (ik | x-x’ |  ) 
JJtt  |  x— x’  | 


(12) 


-  25  - 


We  also  define  the  two-point  correlation  function 


<M(x)y(x'  )>  =  R(x  ,x'  ) 


(13) 


which  is  assumed  to  depend  only  on  x-x*.  With  the  additional 
assumption  of  isotropic  turbulence,  namely  R  =  R(|x-x'|),  the 
Fourier  transform  of  equation  (11)  may  be  written  in  the  form 


0  =  (K2-k2)<|>(K) 


+  k 


G(o)(  |x-x’  |  )R(  |x-x*  |  )<Hx'  )  e1*’*  d3xd3x'  (14) 


where 


4>(K)  = 


it'x  .3 

4> ( x )  e  dJx. 


Finally,  it  is  assumed  that  the  correlation  length  is  small 
compared  to  the  dimensions  of  the  medium.  Then  after  making 
the  transformation,  y  =  x-x’,  the  integrals  of  equation  (14) 
are  separable,  and  the  resulting  dispersion  relation  is  found 
to  be 


0  "  K2  -  k2  + 


exp(ik|y |+iK.y)  +  3 

- - -  H  ( | y  I  )d  y 


(15) 


- - — UntWIftH  I..  II. 


* 
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where  Green’s  function  (12)  has  been  used.  Although  equa¬ 
tion  (15)  has  been  derived  by  assuming  that  k  is  constant 
throughout  the  wake,  it  is  reasonable  to  generalize  this 
equation,  for  a  spatially  dependent  k,  so  that 


0  = 


KOE)?-k(J)2 


+  k(x) 
4tt 


exp[ik(x) |y |+iK(x)*y]  .-►.  3 

- — -  R( I y I )d  y 


y 


(16) 

under  the  condition  that  k(x)  does  not  change  significantly 
in  a  distance  1/k.  The  final  version  of  the  full  model  for 
the  coherent  field,  as  proposed  at  the  Wcrkshop,  involves  the 
replacement  of  equation  (11)  by 


[v2  +  K  ( x ) 2]  <j>  =  0  (17) 

where  K(x)  is  the  solution  of  the  integral  equation  (16). 

Note  that  K(x)  depends  upon  the  mean  electron  density  through 
the  quantity  k(x)  and  upon  the  statistics  of  the  electron 
density  fluctuations  through  the  correlation  function 

R( | x-x ’ | ) . 

The  infinite  space  Green’s  function  of  equation 
(12)  appears  in  equation  (16)  as  the  factor  exp ( iky )/4^y . 

An  obvious  generalization  would  involve  the  replacement  of 


this  factor,  In  equation  (16),  by  a  Green's  function  appro¬ 
priate  to  a  bounded  medium.  This  proposal  was  considered 
at  the  Workshop,  but  It  was  felt  that  such  a  calculation 
would  be  very  difficult  and  should  not  be  attempted  now. 

Equations  ( 1 6 )  and  (17)  constitute  an  extreme 
simplification  of  the  stochastic  problem,  and  it  is  important 
to  study  the  conditions  of  validity.  In  passing  from  equation 

( *0  to  equation  (5)  one  omits  the  term  1^6^  -  which 

2  2  2 
is  of  order  k  y6^,  while  retaining  the  terms  (V  +k  )6\p  and 

2  2 

k  y<j).  The  operator  V  ,  when  applied  to  6\p,  is  expected  to 

2 

be  of  order  1/L  ,  where  L  is  the  correlation  length,  so  that 
the  condition 


2  2 

yk  L  <<  1 

implies  that  the  omitted  term  is  indeed  small  compared  to 
2 

V  6ip.  There  is,  however,  no  a  priori  reason  to  believe  that 

the  term  6 -  <1^6^  =  k^[y<5i|>  -  <(y6ij>^],  which  is  omitted, 

2 

is  much  smaller  than  the  term  k  y<j>,  other  than  the  fact  that 
the  full  model  DWBA-2c  is  an  improvement  over  simpler  models 
when  applied  to  nuclear  scattering  problems.  The  derivation 
of  equation  (16)  from  equation  (11)  regarded  as  heuristic, 
and  no  corresponding  conditions  were  proposed.  The  condition 
that  the  correlation  length  be  small  compared  to  the  dimen¬ 
sions  of  the  wake  is  obvious  . 
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The  "expanded  model"  (DWBA-2b)  is  a  minor 
specialization  of  the  full  model.  The  assumption  is  that 
the  solution  of  equation  (16)  for  K  is  not  very  different 
from  k,  so  that  an  expansion  is  possible.  The  solution  of 
equation  (16)  is  approximately  by 


K  = 


R(y)d3y. 

*y 


(19) 


This  value  of  K  may  be  no  worse  than  the  value  obtained 
from  the  full  model.  The  validity  conditions  are  clearly 
the  same,  since  it  has  already  been  assumed  above  that  the 
fluctuations,  by  means  of  which  K  differs  from  k,  may  be 
treated  as  a  perturbation. 

The  heuristic  model  is  somewhat  of  an  incon¬ 
sistency,  but  was  introduced  in  order  to  obtain  answers  very 
quickly.  Essentially  it  is  the  Born  approximation,  but  with 
the  beam  allowed  to  attenuate  exponentially  with  a  damping 
coefficient  given  by  the  imaginary  part  of  expression  (19), 
that  is,  one  writes  for  the  field 


where 


E  =  Eq  exp(ikx-ax) 


a 


00 

y  sin2(ky)R(y )dy. 


0 
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This  expression  for  a  is  obtained  from  equation  (19)  after 
some  elementary  angular  integrations.  The  heuristic  model 
includes  scattering  out  of  the  beam  to  lowest  order  in  the 
fluctuations,  but  then  it  requires  the  assumption  that  this 
same  value  of  a  can  be  used  to  all  orders  in  an  exponential 
attenuation.  It  ignores  scattering  into  the  beam  (and  may 
therefore,  be  applicable,  in  some  sense,  to  a  thin  beam), 
as  well  as  modifications  due  to  fluctuations  of  the  real 
part  of  the  index  of  refraction. 

The  Kraichnan  model  goes  beyond  the  full  model, 
and  is  the  most  general  perturbation  technique  which  was 
considered.  According  to  the  full  model  of  equation  (7), 
one  can  write  a  generalized  averaged  Green's  function  for 
the  problem  in  the  form 

<0>  -  G(o)  +  G(o)<L1G(°)L1>G(o)  (20) 

whereby  the  statistics  are  included  directly  in  the  Green's 
function.  According  to  the  Kraichnan  model,  equation  (20) 
is  replaced  by  the  nonlinear  integral  equation, 

<a>  -  a(o)  +  o(o)<l1<q>l1Xo> 


(21) 
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No  discussion  was  given  as  to  the  manner  in  which  equation 
(21'>  is  an  improvement  over  equation  (20)  and  possibly  the 
answer  is  unknown.  At  the  present  time,  however,  this  may 
be  an  academic  question,  unless  a  relatively  simple  method 
of  solving  equation  (21)  is  proposed. 

In  all  of  these  methods  the  result  is  a  solution 
for  the  coherent  field  <}>  •  For  a  given  realization  of  the 
random  medium,  that  is,  for  a  given  choice  of  the  function 
=  k(x)  y(x),  the  corresponding  incoherent  field  would  in 
principle  be  calculated  from  equation  (6) 


6tp(x)  =  -Lo^L^4) 


d3x 


’G(o)(x,x’ 


)k(x ’  )2y  (x* ) 4> ( x  T  ) 


(22) 


The  scattered 

<|  6<Kx)  |  2>  = 


energy , 

* 

d3x  * d 


however,  is  proportional  to  the  quantity 
3x"k(x’  )2k(x")2G(o)(x,x’  )G(o)(x,x")* 


•<y(x’  )y(x" ) ( x*  )<J>(x"  )* 

(23) 

2 

Given  a  coherent  field  4>,  the  quantity  <(|<$^|  can  be  calcu¬ 
lated  since  the  correlation  function  <(y  ( x  ’  )y  ( x"  )^>  =  R (  |  x  1  -x"  |  ) 


is  assumed  known. 
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Section  3-C  Watson* s  Transport  Equation 

K.  M.  Watson  has  studied  multiple  scattering  of 
electromagnetic  waves  In  a  random  medium,  starting  from  the 
problem  of  scattering  by  Individual  electrons.  Scattering 
from  a  collection  of  N  scattering  centers  (the  Individual 
electrons)  can  be  exactly  described  by  a  collection  of 
coupled  multiple  scattering  equations.  When  appropriate 
statistical  averages  are  taken  over  the  positions  of  Indi¬ 
vidual  electrons.  It  becomes  possible  to  separate  the 
scattered  waves  into  coherent  and  incoherent  parts.  The 
coherent  wave  may  be  assumed  to  propagate  in  a  medium  of 
varying  effective  wave  number  which,  as  in  the  perturbation 
theories  discussed  earlier,  is  expressed  in  terms  of  the 
electron  density  and  the  correlation  function  of  the  electron 
density  fluctuations. 

Under  suitable  conditions,  the  most  important  being 
the  validity  of  the  eikonal  approximation,  Watson  derives  a 
transport  equation  for  the  intensity  I^j  of  the  incoherently 
scattered  wave: 


dlj_i  (x,0) 

ds 


(ii |M(p ,0 ' ) | tt)dfi^ 


Ijj  (x,p) 


+ 


* 

(ij | M(p ,£  * ) |st)lst(x,e* Jdflg, 


3 


(2*0 


i  =  1,2,  j  =  1,2 
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Here  Iij(x,p)  is  the  four-component  intensity  function, 
corresponding  to  polarization  directions  i , j  =  1,2,  for  the 
incoherent  waves  propagating  in  direction  p  at  position  x. 
The  scattering  function  is  a  4*4  matrix,  which  is  given  by 

(ij  |M(0,0» )  |st)  -  ^  otp2(x)  (s)]  [s~(j  )*^t  (t)] 

in,  (  z)k(p 

•  g(z,R)  e  1  dJR  (25) 

where  aT  is  the  Thompson  scattering  cross  section,  p(x)  is 
the  local  electron  density,  g(z,R)  =  g(z,|z'-z  |)  is  the 
pair  correlation  function  (assumed  to  be  spherically 
symmetric),  and  the  polarization  directions  are  defined  by: 


If  a  solution  to  equation  (24)  for  1^  is  obtained,  then  the 

“►  A 

energy  density  at  a  point  z  with  polarization  e  is  given  by 
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U  ( z ;  6 ) 


1 


<  |  (z)  | ,<L  > 


2 


E  Hp(i)][® 

i,j=l 


(26) 


It  is  clear  from  inspection  of  equation  (22)  that 
the  eikonal  approximation  has  been  made  in  its  derivation, 
except  that  the  wave  properties  associated  with  polarization 
have  been  retained  (and  therefore  1^  has  four  components). 

The  principle  value  of  Watson’s  derivation  lies  in  the  explicit 
evaluation  of  the  scattering  matrix,  equation  (25),  in  terms 
of  the  pair  correlation  function.  In  the  derivation  of  this 
expression  it  is  necessary  to  make  the  same  assumptions  re¬ 
quired  in  the  distorted  wave  Born  approximations. 

According  to  equation  (24)  the  scattering  length 
£(x)  is  given  by 


-V  =  \  S  (ii  |M(£,p'  )  |tt)df2^, 
£(x)  i,t  ■ 


(27) 


=  9  ( x ) 


aT(p  *P* )g(z,R)  e 


in1(x)(^’-p) 


d3RdP 


A  | 

P' 
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where  a,p(£*£')  is  the  differential  Thompson  scattering  cross- 
section.  In  order  of  magnitude,  one  has 

r  *  roRcSp2>  (28) 

2  2 

where  rQ  *  (e  /me  ),  the  classical  electron  radius,  F?c  is 
the  correlation  length,  and  6p  is  the  electron  density 
fluctuation . 

The  most  important  assumptions  required  in  Watson's 
analysis  are 

(i)  Scattering  occurs  in  the  wave  zone,  namely  i  >>  X, 
where  X  is  the  radiation  wave  length; 

(ii)  Only  one  scatter  occurs  in  each  correlated  volume, 
namely  l  >>  Rc ; 

(iii)  Geometrical  optics,  namely  X|V(£n  n)|  <<  1  where 
n  is  the  index  of  refraction; 

(iv)  All  assumptions  required  by  the  distorted  wave  Born 
approximations,  that  is,  sufficiently  underdense 
electron  density  fluctuations; 

(v)  In  order  to  avoid  bending  of  the  ray  paths  (note 
that  s  of  equation  (22)  is  measured  along  a  ray 
path)  as  well  as  to  avoid  rotation  of  polarization 
during  propagation  one  requires  |n-l|  <<  1,  namely 
an  underdense  plasma. 
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Section  3.D  Information  Theory  Formulation 

The  informational  entropy  maximization  principle 
represents  an  attempt  to  obtain  closure  of  the  stochastically 
nonlinear  random  equation  of  wave  propagation,  without  either 
using  a  perturbation  expansion  or  invoking  a  priori  closure 
assumptions.  It  provides  closure  at  any  specified  level  of 
the  hierarchy  of  moment  equations,  by  requiring  that  the 
included  unknown  moments  lead  to  a  maximum  informational 
entropy,  and  therefore  provides  an  unbiased  estimate  of  the 
included  moments  of  the  field,  the  higher  moments  remaining 
unspecified . 

Consider  the  wave  equation 

V2\p  +  k2  (  1+y  )ip  =  0  (29) 

where  p  is  stochastic,  and  not  necessarily  small.  The  first 
two  equations  of  the  moment  hierarchy  are 

(V2+k2)<^(x),>  +  k2<jj  (x)i|/(x )}  =  0 

(v2+k2)<y(x' )^(x)>  +  k2<p(x)p(x' )<Kx)>  =  0  (30) 

The  last  terms  of  both  of  the  equations  (30)  are  the  stochas¬ 
tically  nonlinear  terms  which  lead  to  the  closure  problem. 
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The  maximum  informational  entropy  principle  proceeds  con¬ 
ceptually  as  follows. 

For  any  given  choice  of  <y(x)y(x '  )ip(x)>,  the  equations 
(30)  are  a  closed  system  and  can  in  principle  be  solved  in 
order  to  obtain,  in  addition  to  <y  (x)y(x'  )^(x)*>  ,  the  ensemble 
averages  <y(x' )i|j(x))>  and  <\p(x))>  .  But  by  definition 


<y(x)y(x'  )i|i(x)> 


(3D 


P(y  (x )  =  y,y(x»)  =  y  '  }P(xp(x)  =  iJj  |  y  ,y  '  }yy  '  ipdydy  '  dijj , 


where  P{y(x)  =  y,y(x')  =  y ' }  is  the  joint  probability  that 

y(x)  and  y(x')  will  assume  the  values  y  and  y'  respectively 

and  P  { ip ( x )  =  ^|y,y*}  is  the  conditional  probability  that 
c 

*Hx)  will  assume  the  value  \j>,  given  that  the  values  of  y(x) 
and  y(x')  are  y  and  y*.  The  other  two  ensemble  averages  of 
equations  (30)  are 


<y(x'  )iMx)>  = 


(32) 


P(y(x)  =  y,y(x')  =  y’}P  (ij/(x)  =  ip  |  y  ,y  '  }y  ’i^dydy  T  dip 


and 
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<<Mx)> 


P{y(x)  =  y , y ( x * )  =  y’}P  {if»(x)  =  4>|y,y  1  }i|>dydy ’dif/ 

V 


(33) 

Now  the  probability  P{y(x)  =  y,y(x')  =  y ’ }  is  a  prescribed 
function  which  is  determined  by  the  stochastic  properties 
of  the  medium.  Therefore  by  specifying  the  quantity 
<y(x)y(x'  )i|»(x))> ,  equation  (31)  is  a  constraint  on  the  con¬ 
ditional  probability  P  (ij/(x)  =  4T| pT, TT*  >  .  Furthermore,  the 
wave  equations  (30)  could  then  be  solved  uniquely  for 
<y(x')^(x))>  and  <i^(x)>,  so  that  equations  (32)  and  (33) 

constitute  two  additional  constraints  on  P  .  Given  these 

c 

three  constraints,  there  exists  a  unique  conditional  proba¬ 
bility  function  P  (ij>(x)  =  ty|y,y'}  which  maximizes  the  corre- 

c 

sponding  informational  entropy, 


/ 


Sc(x,x») 


P  _(4»(x)  =  ^|y,y’}£n  P  { ( x )  =  |  y  ,y  '  }di|»  (3*0 

V  O 


The  stationary  value  of  Sc,  say  S  (x,xf),  is  determined  by 
the  condition 


(35) 
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The  resulting  value  of  §c  depends  on  the  original  choice  of 
<Kx,x«)  =  <u(x)y(  x'  )i(/(x));  in  fact  §c  is  a  functional  of  <f>. 
Therefore  maximization  of  §c  with  respect  to  <j> » 

«§  m 

— t*—  -  0.  (36) 

provides  one  with  the  best  choice  of  <j>(x,x’),  consistent  with 
the  retention  of  two  moment  equations  (30).  Clearly  the 
method  can  be  extended  to  an  increased  number  of  moment  equa¬ 
tions,  with  a  great  deal  of  increased  difficulty. 

In  summary,  for  a  given  choice  of  <J>(x,x*)  = 

(x  )y  (x  '  (x  ))>,  equations  (30)  are  solved  for  <m(x'  )\p(x)) 
and<fy(x))>.  Then  S  of  equation  (3*0  is  maximized  with 
respect  to  p  ,  subject  to  the  constraints  (31),  (32),  and 

V 

(33).  Finally,  the  result  is  maximized  with  respect  to  <f> . 

An  increased  number  of  moment  equations  improves  the 
accuracy  by  providing  more  constraints  on  Pc 

The  procedure  described  above,  while  well  defined, 
has  not  been  significantly  explored  at  this  time.  It  is 
discussed  here  because  it  appears  to  have  potential  as  a 
systematic  and  practical  variational  procedure,  and  because 
it  has  been  found  to  be  useful  in  studies  of  equilibrium  and 
nonequilibrium  statistical  mechanics. 


-  39  - 


Section  3.E  Computer  Experiments 

Hochstim  has  performed  a  computer  "experiment" 
which  Is  Intended  to  simulate  scattering  from  a  one¬ 
dimensional  random  medium.  He  solves  the  problem  of  propa¬ 
gation  through  a  stack  of  a  finite  number  of  slabs,  each  of 
which  may  be  assigned  a  thickness  and  a  dielectric  constant. 
The  transmitted  and  reflected  intensities  are  calculated  for 
one  given  set  of  slabs;  this  is  a  well-defined  mathematical 
problem  which  reduces  to  matrix  inversion.  The  calculation 
is  then  repeated  many  times  with  the  thickness  and/or  the 
dielectric  constant  of  each  slab  varying  in  a  random  fashion 
from  one  realization  to  the  next.  These  random  variations 
can  be  described  statistically,  and  the  statistical  properties 
of  an  ensemble  of  such  realizations  are  hence  known.  The 
outputs  of  such  an  "experiment"  for  one  ensemble  of  dielectric 
constant  profiles  are  the  ensemble  averages  of  the  transmitted 
and  reflected  intensities. 

The  "experimental"  results  can  be  compared  with 
theoretical  calculations  for  propagation  in  a  one-dimensional 
random  medium.  The  theoretical  results  consist  of  approxi¬ 
mate  solutions  to  the  one-dimensional  wave  equation  for  a 
medium  with  a  randomly  varying  dielectric  constant;  these 
solutions  can  be  derived  using  various  techniques  which  are 
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analogous  to  the  approaches  to  the  three-dimensional  scat¬ 
tering  problem  (the  first  Born  approximation,  the  distorted 
wave  Born  approximations,  etc.).  Comparisons  of  the  one¬ 
dimensional  theoretical  and  "experimental"  results  give  some 
indication  of  the  ranges  of  validity  of  the  various  theoret¬ 
ical  approaches. 

It  was  proposed  that  there  might  be  several  types 
of  extensions  to  Hochstim's  approach. 

Extensions  of  the  one-dimensional  model.  A  number 
of  additional  features  of  scattering  phenomena  might  be  studied 
within  the  context  of  Hochstim’s  recent  one-dimensional  studies 

a.  Thus  far  Hochstim  has  assumed  a  spatially  uniform 
mean  electron  density  upon  which  t  atistlca.'1  fluc¬ 
tuations  are  superimposed.  Similar  calculations 
could  be  performed  using  a  more  realistic  spatially 
non-uniform  (possibly  parabolic)  mean  electron 
density  distribution. 

b.  Average  phase  shifts  for  the  reflected  and  trans¬ 
mitted  waves  could  be  calculated. 

c.  Time  correlations  between  the  various  statistical 
realization  could  be  specified. 

d.  The  statistics  of  the  random  medium  could  be  made 
more  realistic.  For  example,  Hochstim  assumes  that 


I 
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the  thickness  of  each  slab  and  Its  electron  density 
are  statistically  Independent,  whereas  these  param¬ 
eters  should  be  correlated. 

e.  Hochstim's  results  for  the  reflected  and  trans¬ 
mitted  Intensities  are  presented  as  a  function  of 

2 

the  parameter  (6n/n)  ,  where  6n  Is  the  RMS  electron 
fluctuation  and  n  Is  the  mean  electron  density. 

The  computer  experiments  might  suggest  meaningful 
alternative  parameters. 

f.  Attenuation  of  the  beam  due  to  side  scattering  could 
be  Included  even  In  the  context  of  a  one-dimensional 
calculation. 

g.  The  incident  wave  could  be  allowed  to  impinge  upon 
the  stack  of  slabs  at  an  angle,  in  order  to  study 
aspect  angle  dependence. 

Solutions  of  the  multiple  scattering  equation  of 
Watson.  If  there  are  N  scatterers,  which  are  in  the  wave 
zone  relative  to  each  other  (scatterer  separations  large 
compared  to  the  wave  length),  then  Maxwell's  equations  are 
exactly  replacable  by  the  system  of  algebraic  equations, 


N  2 

!(!„>.  *I«V  £  £  Sa6(J)FaB(5cij) 

8*a«l  j-l 


V 


-  42  - 


and 

G°0fii(ue,eo)Ei(za) 

(37 

N  2 

+  E  E  G«8fiJ(c,8>80)F60<J6J) 

o/B=l  j=l 

'  -+  ^  A 

where  E,(za)  Is  the  Incident  field  at  position  za>eag(J) 
describes  the  polarization  of  the  field  while  travelling 
from  points  z^  to  z^  ;  ^(08,80)  =  _r0^a£  ( 1 )  *  £ga  U  )  ls  the 
Thompson  scattering  amplitude;  f11(aB,80)  Is  the  scattering 
amplitude  for  the  incident  field;  and 


F  0(z  ,1)  = 
a8  a* 


lk I z  -  z  I 
,  1  a  e1 


z  -  z„ 
a  8 


Is  the  free-space  Green's  function  in  the  wave  zone.  The 
solution  of  equations  (37)  involves  diagonallzatlon  of  an 
NxN  matrix.  This  might  be  done  fairly  easily  with  N  as  large 
as  100,  in  order  to  serve  as  a  check  on  scattering  theories, 
and  if  computer  costs  permit,  it  might  be  possible  to  include 
a  statistical  distribution  of  scatterer  realizations. 

3  9 


Laboratory  experiments  on  a  simple  physical  model . 

A  statistical  distribution  of  metal  spheres  might  be  con¬ 
structed,  and  either  microwaves  or  acoustic  waves  scattered 
from  the  distribution.  With  the  statistics  completely  pre¬ 
scribed,  this  is  an  additional  check  on  scattering  theories. 
There  was  some  disagreement  among  members  of  the  Workshop 
concerning  the  feasibility  of  this  experiment. 

One-dimensional  doppler  spread.  Hochstim's  one¬ 
dimensional  calculations  could  be  redone  with  the  boundaries 
of  the  slabs  being  given  random  velocities.  Because  of 
multiple  scattering,  the  transmitted  and  reflected  waves 
would  reflect  a  sequence  of  Doppler  shifts  from  the  individual 
boundaries.  A  question  to  be  answered  by  such  a  study  would 
be  the  effect  upon  the  frequency  spectrum  of  significant 
multiple  scattering. 
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Chapter  4 

LABORATORY  FACILITIES  FOR 
SCATTERING  EXPERIMENTS 

There  have  been  several  laboratory  investigations 
of  the  scattering  of  microwaves  from  turbulent  plasmas.  The 
scattering  media  in  these  experiments  have  been  turbulent 
plasma  Jets  and  ionized  turbulent  pipe  flows.  These  facil¬ 
ities  have  made  possible  direct  comparisons  with  scattering 
theories.  Published  results  of  the  BTL  (Granatstein)  pipe 
flow  experiment  and  the  SRI  (Guthart  et  a] '  plasma  Jet  are 
included  among  the  references  listed  in  Appendix  A. 

In  the  SRI  and  BTL  experiments  the  mean  electron 
density  can  be  varied  from  a  level  of  at  least  two  orders 
of  magnitude  below  critical  electron  density,  up  to  at  least 
critical  electron  density.  The  ratio  of  electron  density 
fluctuations  to  mean  density  is  about  0.5  Tor  both  experiments. 
In  both  facilities  it  is  possible  to  study  the  cross- 
polarized  return  and  the  frequency  spectrum  of  the  scattered 
fields . 

Recently  the  RCA  (Montreal)  laboratory  (Johnston 
et  al)  has  begun  to  study  microwave  scattering  from  a  turbu¬ 
lent  plasma  Jet.  In  their  present  configuration  the  scatter¬ 
ing  is  from  regions  in  which  the  electron  density  exceeds  the 
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critical  electron  density.  (More  recent  results,  obtained 
since  this  Workshop,  have  given  information  about  scattering 
at  lower  electron  densities.) 

In  all  of  these  facilities  there  has  been  a  careful 
diagnostic  study  of  the  properties  of  the  turbulent  medium. 

In  addition  an  extensive  diagnostic  study  of  a  plasma  jet  has 
been  performed  by  Demetriades  at  the  Philco  Corporation. 

The  properties  of  the  scattering  media  and  the  re¬ 
sults  of  the  scattering  experiments  are  summarized  in  the 
following  table. 


Diagnostics  and  Scattering  from  Laboratory  Plasmas 
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Scattering  is  apparently  from  overdense  regions. 


Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


Section  5-A  Conclusions  and  Recommendations  of  the  Theoret¬ 
ical  Panel 

Scattering  theory.  (Small  electron  density  fluc¬ 
tuations.)  In  Chapter  3  we  reviewed  a  number  of  scattering 
theories  which  are  currently  felt  to  be  fruitful.  In  the 
limit  of  small  electron  density  fluctuations  the  first  Born 
approximation  is  frequently  applicable,  as  well  as  extensions 
to  it  such  as  the  various  distorted  wave  Born  approximations. 
The  panel  recommends  that  extensions  of  the  first  Born  approx¬ 
imation  should  continue  to  receive  extensive  study,  especially 
with  regard  to  applications  to  the  wake  scattering  problem; 
it  is  felt  that  the  "full  model"  (DWBA-2c)  should  receive  the 
greatest  attention.  The  full  model  is  chosen  over  some  of 
the  simpler  extensions  discussed  in  Chapter  3  (namely  DWBA-1, 
-2a,  and  -2b)  because  it  is  only  slightly  more  complicated 
and  is  expected  to  yield  significantly  better  results. 

It  appears  (on  the  basis  of  the  studies  presented  at 
the  Workshop)  that  the  first  Born  approximation  has  a  greater 
range  of  validity  than  had  been  previously  acknowledged.  The 
limits  of  validity  of  the  Born  approximation,  including  the 
Salpeter-Treiman  condition,  should  be  reconsidered. 
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The  role  of  Watson’s  transport  theory  as  a  supple¬ 
ment  to  the  various  perturbation  approximations  should  be 
studied  further. 

Scattering  theory.  (Partially  or  totally  overdense 
media.)  Studies  of  scattering  from  overdense  random  surfaces 

I 

should  be  continued  and  the  theories  modified,  if  necessary, 
in  order  to  make  them  applicable  to  scattering  from  an  over- 
dense  turbulent  wake. 

As  a  preliminary  to  the  study  of  scattering  from  a 
wake  with  overdense  blobs,  scattering  and  absorption  by 
single  overdense  geometrical  structures  should  be  studied.  A 
considerable  literature  on  this  subject  is  available  in  the 
context  of  ionospheric  physics.  It  is  known  that  considerable 
absorption  occurs  when  an  electron  density  profile  becomes 
overdense  with  a  gentle  slope;  therefore,  it  is  conceivable 
that  the  inclusion  of  absorption  may  simplify  the  multiple 
scattering  problem  by  reducing  the  number  of  scatterers  re¬ 
quired  for  the  investigation. 

Also  a  preliminary  to  the  study  of  overdense 
scattering,  it  may  be  useful  to  develop  techniques  for  ray 
tracing  with  representative  electron  density  profiles. 

The  panel  is  not  optimistic  about  the  chances  of 
success  for  these  approaches  to  scattering  from  partially  or 
totally  overdense  media. 
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Inputs  to  theory.  (Representation  of  the  scatter¬ 
ing  medium.)  It  is  very  important  to  develop  a  mathematical 
description  of  a  random  scattering  medium  which  is  clearer 
than  is  presently  available  in  the  statistical  descriptions. 
Thus  if  a  wake  could  be  characterized  as  an  irregular  tube 
with  sharply  defined  boundaries  within  which  the  electron 
density  is  reasonably  uniform  and  outside  of  which  the  elec¬ 
tron  density  is  zero  the  proper  electromagnetic  theory  to 
use  would  most  likely  be  one  in  which  surface  scattering  was 
equally  as  important  as  volume  scattering.  On  the  other 
hand,  the  statistical  description  of  electron  density  for  a 
wake  could  imply  a  theory  such  as  a  distorted  wave  Born  approx¬ 
imation.  Therefore,  it  is  recommended  that  fluid  dynamicists 
make  an  effort  to  provide  a  description  of  a  single  realiza¬ 
tion  of  a  turbulent  wake  in  sufficient  detail  to  enable 
scattering  theorists  to  choose  the  correct  model  of  scattering. 

Inputs  to  theory.  (Scattering  experiments.)  A 
major  problem  in  the  development  of  scattering  theory  is  the 
determination  of  the  effective  index  of  refraction  (or  equiv¬ 
alently,  the  effective  wave  number)  in  a  turbulent  medium. 

The  Importance  of  this  information  is  clear  from  the  review 
of  scattering  theories  in  Chapter  3.  It  is  therefore  recom¬ 
mended  that  microwave  scattering  experiments  be  designed  to 
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measure  the  forward  scattered  wave  ,  as  a  function  of  scatter¬ 
ing  angle,  as  well  as  the  associated  phase  shifts,  and  that 
these  experiments  be  analyzed  to  determine  the  propagation 
properties  of  the  turbulent  medium. 

Studies  of  frequency  spectra.  Studies  of  the 
frequency  spectra  of  electromagnetic  waves  scattered  from 
turbulent  media  have  appeared  in  the  literature,  almost 
exclusively  in  the  context  of  the  first  Born  approximation. 
However,  the  effect  of  multiple  scattering  upon  the  frequency 
spectrum  is  an  extremely  important  problem  upon  which  there 
is  presently  little  theoretical  information.  Therefore,  the 
current  effort  should  be  extended  to  calculations  of  fre¬ 
quency  spectra  using  the  distorted  wave  Born  approximations. 

Numerical  experiments.  The  feasibility  of  extend¬ 
ing  Hochstim's  computer  experiment  to  three  dimensions  should 
be  investigated.  For  each  realization  of  an  ensemble  of  ran¬ 
dom  distributions  of  scatterers  this  "experiment"  would 
involve  the  diagonalization  of  an  N*N  matrix,  where  N  is  the 
number  of  scatterers  in  the  system.  The  study  should  begin 
by  calculating  the  scattering  from  a  single  realization, 
thereby  modeling  the  scattering  problem  without  statistics. 

If  the  cost  is  reasonable,  then  many  realizations,  perhaps 
with  specified  statistics  (or  at  least  known  statistics), 


should  be  studied.  It  is  recommended  that  research  programs 
along  these  lines  have  relatively  low  priority. 


Other  useful  extensions  of  Hoch3tlm's  work,  all 
in  the  context  of  the  one-dimens lonal  problem,  have  been 
discussed  in  Section  D  of  Chapter  3. 

Some  consideration  should  be  given  to  the  validity  of 
extrapolating  to  three  dimensions  the  results  of  one-  and 
two-dimensional  numerical  experiments. 
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Section  5*B  Conclusions  and  Recommendations  of  the  Experi¬ 
mental  Panel 

Scattering  mechanisms.  The  first  Born  approximation 
for  scattering  from  a  turbulent  plasma  appears  co  be  valid  In 
the  far  wake  of  many  of  the  ballistic  range  projectiles  and 
full-scale  reentry  vehicles  which  have  been  studied.  The 
range  of  applicability  of  the  first  Born  approximation  appears 
to  be  greater  than  that  given  by  the  Salpeter-Treiman 
conditions . 

Departures  from  the  predictions  of  the  first  Born 
approximation  may  be  caused  by  one  or  more  of  the  following 
physical  processes: 

a)  The  .presence  of  locally  overdense  "blobs"  in  an 
otherwise  underdense  plasma. 

b)  Attenuation  due  to  absorption  and/or  scattering. 

c)  Variation  of  the  effective  wave  number  of  the 
coherent  wave  propagating  in  the  turbulent  medium. 

The  mechanism  (a)  was  not  considered  in  any  detail  by  the 
Workshop;  mechanisms  (b)  and  (c)  are  taken  into  account,  at 
least  to  first  order  in  the  electron  density  fluctuations 
in  the  various  distorted  wave  Born  approximations.  The  rela¬ 
tive  importance  of  these  various  mechanisms  has  not  been 
specifically  studied.  Only  volume  scattering  effects  were 
considered  by  the  Workshop;  no  attention  was  given  to  surface 
scattering. 
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Comparisons  of  theory  and  experimental  -data.  For 
sufficiently  small  electron  density  fluctuations  scattering 
from  turbulent  laboratory  plasmas  appears  to  agree  with  the 
first  Born  approximation,  departures  from  the  first  Born 
approximation  (such  as  those  seen  in  the  experiments  of 
Granatstein)  may  be  satisfactorily  interpreted  in  terms  of 
attenuation  effects,  as  in  the  theories  of  Shkarofsky  and  of 
Feinstein.  (These  theories  fall  into  the  category  DWBA-2a 
of  Chapter  3 .  ) 

Preliminary  attempts  were  made  (during  the  Workshop) 
to  interpret  a  limited  amount  of  wake  scattering  field  data 
from  non-ablating,  non-seeded  spheres  at  high  altitudes.  It 
appears  possible  to  achieve  a  self-consistent  interpretation 
of  this  data  using  (i)  the  first  Born  approximation,  (ii) 
attenuation  corrections  to  the  first  Born  approximation  using 
the  heuristic  models  (DWBA-2a)  of  Shkarofsky  and  of  Feinstein, 
(iii)  extrapolation  to  reentry  conditions  to  turbulence 
theories  and  experiments  on  low  speed  wakes  and  plasma  jets, 
assuming  local  isotropy  and  uniform  molecular  transport  prop¬ 
erties.  These  interpretations  must  be  considered  preliminary 
until  more  reliable  values  of  electron  densities,  electron 
density  fluctuations,  and  turbulent  scale  lengths  become 


available . 
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Similar  attempts  were  made  to  interpret  wake 
scattering  field  data  for  heavily  ablating  reentry  vehicles. 
These  attempts  were  not  successful.  It  appears  that  for  many 
ablating  reentry  vehicles  the  electron  densities  are  suffi¬ 
ciently  large  that  strong  attentuation  and/or  multiple 
scattering  effects  must  be  considered. 

Studies  of  this  sort  should  be  continued,  with  an 
effort  being  made  to  use  the  best  available  estimates  of 
turbulent  spectra.  Using  presently  available  theoretical 
models  (one  of  the  distorted  wave  Born  approximations,  for 
example)  field  data  should  be  tested  for  a  consistent 
interpretation  in  terms  of  these  models. 

Studies  of  the  nature  of  the  wake.  The  attempts 
to  obtain  a  consistent  interpretation  of  field  data  have 
suggested  areas  where  further  theoretical  and  experimental 
studies  of  the  properties  of  wake  turbulence  are  necessary. 

a)  For  any  given  reentry  vehicle  there  are  significant 
variations  in  the  turbulence  spectrum  (including 
the  turbulent  scale  sizes)  as  a  function  of 
Reynolds  number  and  of  downstream  distance.  More 
information  about  these  variations  should  be  ob¬ 
tained  from  theory  and  experiment. 

b)  It  has  been  assumed  during  the  course  of  the 
studies  of  the  field  data  that  the  spectrum  of  the 
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electron  density  fluctuations  Is  Identical  to  the 
spectrum  of  the  fluid  velocity  fluctuations.  In 
the  near  future  this  assumption  will  be  tested 
experimentally  for  various  plasma  conditions, 
c)  It  has  been  assumed  that  the  spectrum  of  electron 
density  fluctuations  Is  Isotropic.  Further  studies 
of  departures  from  Isotropy  In  reentry  vehicle 
wakes  Is  required. 

Many  of  these  studies  of  wake  turbulence  can  be  furthered  by 
laboratory  experiments  In  plasma  jets  and  In  ballistic  ranges 
and  shock  tubes.  Specific  proposals  for  useful  studies  In 
such  facilities  are  given  below. 

Laboratory  experiments.  To  gain  better  under¬ 
standing  of  the  properties  of  turbulent  shear  flows  In  a 
compressible  fluid  with  variable  molecular  transport  co¬ 
efficients,  finite  reaction  rates,  etc.,  experiments  should 
be  designed  to  determine  such  relevant  quantities  as: 

a)  Lateral  and  transverse  correlation  functions  (that 
Is,  the  full  three-dimensional  wave  number  spectrum) 
of  velocity,  temperature,  and  electron  density 
fluctuations  and  the  dependence  of  these  correlation 
functions  on  all  of  the  various  parameters  of 
Interest:  Reynolds  number,  temperature  fluctuation 
amplitude,  Debye  length  to  wavelength  ratio,  etc. 
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b)  Spatial  distributions  of  mean  electron  density 
and  of  electron  density  fluctuation  amplitudes  in 
turbulent  shear  flows,  measured  either  directly 
or  through  inference  from  distributions  of  other 
convected  scalar  quantities. 

c)  The  relationship  between  the  spectra  of  electron 
density  fluctuations  and  other  scalar  quantities 
and  the  velocity  spectrum,  over  the  range  of 

parameters  of  interest. 

I 

To  aid  the  study  and  testing  of  scattering  theories 
the  range  of  parameters  under  which  scattering  experiments 
are  performed  should  be  extended.  In  all  of  these  scattering 
experiments  it  is  important  to  adequately  map  the  turbulence 
field  in  the  scattering  medium. 

All  of  these  laboratory  investigations  should  be 
attempted  in  ballistic  ranges  or  shock  tunnels  so  that  the 
differences  between  low  velocity  and  hypersonically  generated 
wakes  can  be  studied. 

Ballistic  range  and  shock  tunnel  experiments  should 
study  the  variations  in  electron  density  fluctuations  for  a 
fixed  flow  field  profile  in  order  to  determine  quantitatively 
whether  departures  from  the  first  Born  approximation  are  to 
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be  expectec .  Scattering  experiments  on  wakes  in  such 
facilities  should  be  performed  to  demonstrate  quantita¬ 
tively  agreement  or  disagreement  with  various  theoretical 
models  for  scattering. 


] 
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1.  INTRODUCTION 


Interpretation  of  radar  returns  from  turbulent 

plasmas  is  often  complicated  by  such  important  effects  as 

1  2 

multiple  scattering  *  and  three  dimensional  geometry. 

However,  under  certain  restrictive  conditions  regarding  the 

plasma  density  and  the  overall  dimension  of  the  scattering 

volume,1  it  is  generally  agreed  that  the  averaged  scattering 

intensity  can  be  reliably  calculated  according  to  the  Booker 

formula,  which  was  a  single-scattering  theory  based  on  the 

4 

first-order  Born  approximation  and  requires  only  relatively 
simple  specification  of  the  statistical  properties  of  the 
scattering  medium.  The  Booker  formula  gives,  for  the  time- 
averaged  scattering  intensity  from  a  turbulent  plasma  of 
sufficiently  slowly  varying  spatial  configuration,  ,:> 


4  2 

./«  \  _  k  sin  \p 

(0>V  — 


Ae^(5;)F(5,g)d3r 


(1) 


V 


Here  a(0,AQ)  is  defined,  as  usual,  as  4tt  times  the  power 
scattered  per  unit  solid  angle  per  unit  incident  power 
density  along  the  direction  of  scattering  which  makes  an 
angle  0  with  respect  to  the  incident  wave  propagation  vector 
^o;k  =  |^|  =  |^o|  =  27T /X ^  is  the  incident  wave  number;  ip  is 
the  angle  between  £  and  the  electric  field  vector  of  the 
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2 

incident  wave;  Ae  (£)  is  the  local  mean-square  fluctuation 
of  the  dielectric  constant  about  its  mean  value  e (5)  at  point 
and 


* 

S(;,5»)  e  dV 
■  V 


(2) 


2 

defines  the  local  spectrum  function  for  Ae  (r)  as  a  spatial 
Fourier  transform  of  the  normalized  two-point  correlation  fun 
tion 

S(£,£')  =  Ae (r )  Ae (r+r ' ) /  Ae2(r)  (3) 

along  the  direction  of  the  vector  <3.  =  £  —  £Q.  The  volume  of 
integration  V  in  equations  (1)  and  (2)  is  understood  to  be 
extended  over  the  effective  scattering  volume  as  defined  by 
the  plasma  boundaries  and/or  the  range  cell  of  the  radar 
under  consideration.  The  frequency-  or  wavelength-dependence 
of  the  scattering  intensity  is,  of  course,  implicitly  con¬ 
tained  in  the  Fourier  parameter  5,  which  has  a  magnitude 
q  =  (1lir/Xo)sin(8/2)  . 

In  the  study  of  radar  return  from  the  turbulent 
wake  of  hypersonic  objects,  there  has  been  a  strong  tempta¬ 
tion  on  the  part  of  some  early  investigators  to  compare  the 
frequency-dependence  of  the  scattering  intensity  o(9,Ao) 
with  what  could  be  Inferred  from  simple  turbulence  theories. 
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In  particular,  the  temptation  to  draw  conclusions  from  a 

limited  number  of  frequency  samplings  about  whether  a 

"universal  spectrum"  of  the  classical  Kolmogorov^  or 

6  7 

Oboukhov-Corrsin  *  form  does  exist  or  not  in  the  turbulent 
hypersonic  wake  seemed  irresistable .  Such  an  exercise  was 
bound  to  be  unfruitful  and  misleading  since  it  not  only  over¬ 
looked  the  fact  that  the  statistical  properties  of  the  tur¬ 
bulent  hypersonic  wake  plasma  within  the  scattering  volume  V 
as  sampled  by  the  finite  range  cell  of  the  radar  might  not 
at  all  be  isotropic  nor  spatially  homogeneous,  but  also  the 
fact  that  even  for  a  homogeneous,  isotropic  turbulence  field 
of  finite  Reynolds  number  in  a  constant-density  fluid  of 
constant  transport  properties  (i.e.,  kinematic  viscosity  and 
molecular  dif fusivity ) ,  a  "universal  spectrum"  can  exist  onlj 
over  a  limited  range  of  wave  number  which,  in  turn,  depends 
on  a  number  of  scaling  parameters . * 10  Furthermore,  the 
effects  of  variable  density  and  rapid  chemical  reactions 
within  the  hypersonic  wake  plasma  may  greatly  complicate  the 
form  of  the  scalar  fluctuation  spectrum . ^  We  shall  review 
these  complications  briefly  as  follows. 

2.  ENERGY  SPECTRUM  FOR  ISOTROPIC  TURBULENCE  IN  AN 

INCOMPRESSIBLE  FLUID  OF  CONSTANT  KINEMATIC  VISCOSITY 

12-] 

As  discussed  at  length  in  most  standard  texts, 
the  three-dimensional  energy  spectrum  E(q,t)  for  a  decaying 
homogeneous  isotropic  turbulence  field  in  an  incompressible 
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fluid  of  constant  kinematic  viscosity  v  is  generally  of  the 
form  illustrated  in  Pig.  1.  As  usual,  the  function  E(q,t) 
is  defined  in  such  a  way  that  the  ensemble-averaged  turbu¬ 
lence  kinetic  energy  per  unit  mass  lying  within  the  wave- 
number  range  between  q  and  q  +  dq  at  time  t  is  given  by 
E(q,t)dq,  so  that  the  mean-square  velocity  fluctuation  along 
any  one  of  the  three  orthogonal  rectilinear  coordinates  is 
given  by  the  integral 


|  u'2(t) 


r 00 

E(q,t)dq  . 

■  0 


(2) 


The  entire  energy  spectrum  is  roughly  divisible  into  three 
main  ranges  In  wavenumber  space.  Referring  to  Fig.  1,  these 
are : 

(i)  The  Low  Wavenumber  Range,  consisting  of  the  largest 

* 

eddies  with  longest  persistency,  and  those  large 
eddies  which  are  mostly  responsible  for  the  macro¬ 
scopic  diffusion  property  (i.e.,  eddy  diffusivity 
Dg)  of  the  turbulence  field.  ^ 


It  should  be  noted  that  "eddies"  in  turbulence  is  actually 
a  loose  term  referring  to  certain  Fourier  components  in  con 
figuration  space.  Eddies  of  different  sizes  are  not  really 
separable  entities  since  they  may  share  the  same  part  of  th 
fluid  at  any  given  time. 
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The  shape  for  this  part  of  the  spectrum  is  slowly 
time-varying,  and  depends  strongly  on  the  initial 
condition  of  formation. 

(ii)  The  Energy-Containing  Wavenumber  Range,  which 

forms  the  hump  of  the  spectrum  in  the  vicinity  of 
the  wavenumber  denoted  by  q£  in  Pig.  1. 

(iii)  The  Universal  Equilibrium  Range  at  high  wavenumber 
beyond  the  hump  of  the  spectrum.  The  shape  of  this 
part  of  the  spectrum  becomes  "universal"  and  inde¬ 
pendent  of  the  initial  condition  of  formation  only 
in  the  sense  that  it  can  be  collapsed  into  a  single 
curve  at  all  time  t  when  the  wavenumber  scale  is 

normalized  with  respect  to  the  Kolmogorov  wave- 
15 

number 


qK  =  (e/v3)*  (3) 

Here  £  denotes  the  instantaneous  rate  of  dissipa¬ 
tion  of  turbulence  energy  per  unit  mass,  which  is, 
of  course,  a  function  of  time  t  in  a  decaying 
turbulence  field. 

It  is  important  to  note  that  the  famous  Kolmogorov  law 
which  predicted  a  q  dependence  for  the  energy  spectrum,  is 
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only  applicable  to  the  "inertial  subrange"  which  lies  some¬ 
where  within  the  lower  wavenumber  portion  of  the  universal 
equilibrium  range  as  illustrated  in  Fig.  1.  In  the  high 
wavenumber  end  near  q^,  the  energy  spectrum  decreases  much 
more  rapidly  than  q  due  to  the  action  of  viscosity.  The 
extent  of  the  inertial  subrange  thus  depends  on  the  instan¬ 
taneous  separation  between  the  two  wavenumbers  q£,  qK  in  a 
decaying  turbulence  field. 

The  inertial  subrange,  as  well  as  the  universal 

equilibrium  range  of  the  energy  spectrum,  has  been  found  to 

12 

exist  not  only  in  low  speed  grid  turbulence,  but  also  in 

IS 

a  large  number  of  shear  flows  ranging  from  tidal  channels 
to  low  speed  wakes  in  water  tunnels. ^  The  normalized 
energy  spectrum  as  deduced  from  these  latter  experiments'*''’*  ^ 

-V, 

showing  the  q  a inertial  subrange  and  the  subsequent  viscous 
cutoff  near  q^  is  reproduced  here  in  Fig.  2. 

As  mentioned  earlier,  the  lower  wavenumber  portion 
of  the  energy  spectrum  (i.e.,  ranges  (i)  and  (ii)  cited  above) 
is  not  universal  even  in  low  speed  grid-generated  turbulence 
of  constant  mean  flow  velocity  U.  This  latter  type  of  flow 
is,  of  course,  a  close  laboratory  simulation  of  a  decaying 
homogeneous  isotropic  turbulence  field  in  which  the  time 
after  formation  t  is  replaced  by  the  averaged  flow  time  past 
the  grid  x/U  at  any  given  axial  station  of  distance  x 
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downstream  of  the  grid.  A  typical  time  evolution  of  the 
normalized  energy  spectrum  in  grid  turbulence  is  illustrated 
in  Pig.  3-  The  lack  of  universality  for  the  lower  wave- 
number  portion  of  the  spectrum  is  quite  evident.  In  the 
case  of  turbulent  shear  flow,  one  may  expect  not  only  simi¬ 
lar  continuous  evolution  of  the  spectral  shape,  but  also 
pronounced  anisotropy  associated  with  the  lower  wavenumber 
portion  of  the  energy  spectrum  as  well. 

3.  RELATIONSHIP  BETWEEN  ENERGY  SPECTRUM 
AND  SCALAR  FLUCTUATION  SPECTRUM 

From  the  Booker  formula  (1),  it  is  quite  clear  that 
the  averaged  scattering  intensity  (or  radar  cross-section) 
does  not  depend  directly  on  the  energy  spectrum  E(£,<j)  of  the 
turbulence  field,  but  rather  on  the  spectrum  function  F(£,g) 
for  mean-square  fluctuation  of  the  dielectric  constant  as 
defined  in  equations  (2)  and  (3).  Even  though  in  subsonic 
turbulence,  the  fluctuation  of  any  scalar  quantity,  such  as 
the  dielectric  constant  e,  is  mostly  driven  by  the  convective 
motion  of  the  turbulence,  there  exists  no  simple  relationship 
between  E(j,g)  and  F(^,g)  in  turbulent  shear  flow. 

In  the  case  of  homogeneous  isotropic  turbulence  of 
sufficiently  high  Reynolds  number  in  a  fluid  of  constant  mass 
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density  p,  kinematic  viscosity  v,  and  molecular  diffusivity 

D  for  a  particular  conserved,  passive  scalar  quantity  0, 

however,  it  has  been  shown  by  Oboukhov,^  and  independently 
7 

by  Corrsin,  that  the  fluctuation  spectrum  for  0  in  the 
"convection  subrange"  should  be  of  the  form 


f(q,t) 


(4) 


in  which  x  denotes  the  instantaneous  rate  of  diffusive 
dissipation  of  mean-square  fluctuation  of  0  per  unit  volume; 
and  6  denotes  the  instantaneous  rate  of  viscous  dissipation 
of  turbulence  kinetic  energy  per  unit  mass,  as  before.  The 
fluctuation  spectrum  T(q,t)  is  defined  in  such  a  way  that 
the  Fourier  component  for  mean-square  fluctuation  of  0  lying 
within  the  wavenumber  range  between  q  and  q  +  dq  at  time  t 
is  given  by  r(q,t)dq,  so  that  the  total  mean-square  (i.e., 
ensemble-averaged)  fluctuation  of  0  at  time  t  is  given  by 


A02(t) 


T(q,t)dq  . 

0 


(5) 


In  a  way  analogous  to  the  definition  of  the  "inertial  sub¬ 
range"  for  the  nergy  spectrum  (Fig.  1),  the  "convection 
subrange"  here  refers  to  a  range  of  wavenumber  lying 
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somewhf  e  between  qg  where  the  peak  of  r(q,t)  is  located, 
and  a  certain  cut-off  wavenumber  which  we  shall  call  the 
Oboukhov-Corrsin  wavenumber, 

qoo  3  (e/D3)*  •  (6 

6  7 

While  Oboukhov  and  Corrsin  both  predicted  that  the  wave- 
number  dependence  of  T(q,t)  within  the  convection  subrange 
should  be  identical  to  that  of  the  energy  spectrum  within 
the  inertial  subrange,  and  that  the  respective  diffusive 
and  viscous  cut-off  wavenumbers  for  the  two  subranges 
should  bear  the  ratio 


(7) 


g 

Batchelor  pointed  out  that  such  would  be  the  case  only  when 
D  and  v  are  of  the  same  order  of  magnitude.  When  D  <<  v, 
such  as  diffusion  of  electrolytes  (e.g.,  salt)  in  water, 

_5/ 

the  q  3  dependence  for  r(q,t)  should  extend  only  to  the 
neighborhood  of  the  Kolmogorov  wavenumber  q^.  In  the  range 
of  wavenumber  between  q^  and  a  new  diffusive  cut-off  wave- 
number  which  we  shall  call  the  Batchelor  wavenumber 
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'5B-(eA®2)  (8) 

the  wavenumber  dependence  of  r(q,t)  should  become  q 

The  existence  of  the  viscous-convective  subrange 

with  q”1  dependence  for  the  fluctuation  spectrum  of  a  weakly 

8 

diffusive  passive  scalar  as  predicted  by  Batchelor  has 
indeed  been  confirmed  in  the  grid  experiment  of  Gibson  and 

"I  *7 

Schwarz,  and  also  more  recently  in  the  wake  experiments 

1  8 

of  Gibson,  Lyons,  and  Hirschsohn.  It  is  worth  noting, 

_5/ 

however,  that  the  transition  point  from  the  q-  3  spectrum  to 
the  q-1  spectrum  observed  in  the  experiments1^* turned  out 
to  lie  in  the  neighborhood  of  0.04  q^  instead  of  q^. 

j| 

When  the  molecular  dlffusivity  for  the  passive 

scalar  of  interest  is  much  greater  than  the  kinematic  vis- 

o 

cosity,  i.e.,  D  >>  v,  Batchelor,  Howells,  and  Townsend"^ 

$ 

predicted  that  the  diffusion  cut-off  to  the  scalar  fluctu¬ 
ation  spectrum  should  begin  at  the  Oboukhov-Corrsin  wave- 
number  qQC,  and  that  in  the  "inertial-diffusive  subrange" 
between  qQC  and  q^,  the  spectrum  should  be  of  the  form 


f(q) 


=  i  Cxe2/3D-V^3 


(9) 


Appendix  B  -  12 


In  other  words,  the  spectrum  should  fall  off  rapidly  with 
Increasing  wavenumber  at  a  rate  approaching 

In  a  series  of  recent  papers,  Gibson‘S  has  pro¬ 
posed  a  unified  spectral  theory  for  the  fine  structure  of 
scalar  fields  mixed  by  turbulence  In  which  the  Inertial  sub¬ 
range  of  Oboukhov  and  Corrsln  «  q~  /s)  and  the  viscous- 
convective  subrange  of  Batchelor  «  q”1^  were  reproduced. 
However,  for  strongly  diffusive  scalars  (i.e.,  D  >>  v), 

Gibson  predicted  a  new  inertial-diffusive  subrange  with 
r  «  q  and  upper  cut-off  at  the  Batchelor  wavenumber  qB» 

This  latter  result  by  Gibson  thus  contradicted  sharply  with 

the  inertial-diffusive  subrange  obtained  earlier  by 

0 

Batchelor:,  Howells,  and  Townsend.  As  of  this  time,  there 
yet  exists  no  definitive  experimental  measurement  on  the 
spectral  structure  of  strongly  diffusive  scalars  in  iso¬ 
tropic  turbulence  from  which  the  shape  of  the  inertial- 
diffusive  subrange  can  be  clearly  established.  In  view  of 
the  fact  that  the  physical  models  employed  by  Batchelor  et 
al.  and  by  Gibson  in  arriving  at  their  respective  theoretical 
results  were  quite  different  ,  and  they  were  believed  to  be 
respectively  correct  with  equal  conviction,  the  question  of 
spectral  shape  for  strongly  diffusive  scalars  in  the  inertial- 
diffusive  subrange  must  be  regarded  as  an  unsettled  problem 
in  basic  turbulent  scalar  mixing  theory. 
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At  this  point,  it  may  be  noted  that  the  scalar 
fluctuation  spectrum  T(q)  discussed  above  is  still  quite 
distinct  from  the  scattering  spectrum  F(g,g)  defined  earlier 
in  Section  1.  However,  it  is  easy  to  show  that  the  two  are 
simply  related  as  follows. 

Assuming  that  the  correlation  length  scale  for 
the  fluctuating  scalar  quantity  of  interest  (e.g.,  the  di¬ 
electric  constant  of  the  wake  plasma)  is  generally  much 
smaller  than  the  dimensions  of  the  scattering  volume  V  under 
consideration  one  may  extend  the  finite  range  of  volume 
integration  in  equation  (2)  to  infinity  without  making 
appreciable  error.  Then,  from  the  well  known  properties  of 
Fourier  transforms,  one  obtains  from  equation  (2), 


s(5',5’) 


8tr ' 


-ig-g'  . 

F(g,g)  e  d Jq 


(10) 


2 

Upon  multiplying  both  sides  of  the  above  equation  Ae  (5) 
and  subsequently  letting  r1  =  0,  one  obtains,  with  the  help 
of  equation  (3)s 

00 

ret  _ 

Ae2(£)R(£,g)d3q  .  (11) 


Ae2(£)  = 

8tt  - 


.00 
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If  the  scalar  field  were  locally  isotropic  about  r,  the 
function  F  varies  only  with  the  magnitude  of  the  vector  ^ 
in  accordance  with  equation  (2),  so  that  equation  (11)  re¬ 
duces  to 


Ae2(r) 


00 

Ae2(r)F(r,q)q2dq  . 

■  0 


(12) 


Comparing  this  with  the  definition  of  the  local  scalar 
fluctuation  spectrum  T  defined  earlier  in  equation  (5), 
with  0  presently  replaced  by  e,  and  t  by  r,  we  obtain  the 
relationship , 


F(r,q)  =  2 it2  T(r,q)/q2  Ae2(r)  (13) 

4.  DIELECTRIC  CONSTANT  AND  ELECTRON  DIFFUSIVITY 
IN  HYPERSONIC  WAKE  PLASMAS 

As  was  discussed  at  considerable  length  in  Ref.  11, 
the  local  dielectric  constant  e(r)  in  a  hypersonic  wake 
plasma  is  generally  given  by  a  linear  combination  of  terms 
involving  the  local  number  density  of  free  electrons  ng(rj, 
the  number  densities  of  molecular  ions  of  both 
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positive  and  negative  charges  n^(r),  and  the  mass  density 

p(r).  Accordingly,  the  two-point  auto-correlation  function 

for  dielectric  constant  fluctuation  S(r.r')  is  generally 

-*■  -► 

made  up  of  a  large  number  of  terms  corresponding  to  auto¬ 
correlation  of,  and  cross-correlation  among,  the  various 
scalar  quantities  just  mentioned.  The  relative  importance 
of  the  various  correlation  terms  in  contributing  to  the 
resultant  dielectric  constant  fluctuation  spectrum  has  been 

shown'*''*'  to  depend  strongly  on  the  electron  chemistry  within 
20 

the  wake,  and  also  on  the  incident  electromagnetic  wave 
frequency  under  consideration. 

In  the  present  discussion,  we  shall  limit  our¬ 
selves  to  the  simpler  type  of  situation  in  which  the  dielec¬ 
tric  constant  fluctuation  is  dominated  by  the  fluctuation  of 
a  single  scalar  quantity,  namely,  the  electron  number 
density  ng.  Such  would  be  the  situation,  for  example,  in  a 
"recombination-controlled’'  wake  where  the  rate  of  disappear¬ 
ance  of  free  electrons  is  governed  by  the  charge-neutraliza¬ 
tion  process  rather  than  by  the  electron-attachment  (negative 
ion  formation)  process.  *  ’  If  one  further  assumes  that 

the  characteristic  time  for  depletion  of  free  electrons  due 
to  chemistry  is  relatively  long  in  comparison  with  the  local 
characteristic  time  scale  for  turbulent  mixing,  then  one  may 
treat  ng  as  a  passive  scalar  quantity  in  the  application  of 
existing  turbulent  mixing  theories. 
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As  is  well  known  in  plasma  physics,  the  effective 

molecular  diffusivity  for  electrons  in  a  weakly-ionized 

plasma  is  a  strong  function  of  the  ratio  between  the  local 
22 

Debye  length 


h  =  (kT/llTTne  e2Y  (1*0 

and  the  local  electron  density  gradient  length  scale  £  of 

interest.  When  h  <<  £,  the  diffusion  process  is  essentially 
21 

ambipolar,  and  the  effective  diffusivity  D  for  the  electron- 
ion  pair  through  the  gas  can  be  expected  to  be  roughly  the 
same  as  the  local  kinematic  viscosity  v  in  a  single  temper¬ 
ature  plasma.  When  h  >>  £,  the  diffusion  of  electrons 
becomes  "free"  from  the  impeding  effect  of  the  heavier  ions, 
and  the  effective  diffusivity  for  these  "free"  electrons 
will  become  much  greater  than  the  local  kinematic  viscosity 
(by  a  factor  of  the  order  of  i.e.,  the  square-root 

of  the  molecule/electron  mass  ratio,  or  a  factor  of  about 
200  in  air).  Thus,  according  to  equation  (1*0,  and  the 
physical  constants  k  =  1.380*4x10”^  erg/°K  (Boltzmann  con¬ 
stant),  e  =  *1 .  803x10""  e.s.u.  (electron  charge),  the  criteria 
for  ambipolar-  and  for  free-dif fusion  of  electrons  are  respec¬ 
tively  given  by  (with  T  in  K  and  £  in  cm). 
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For  amblpolar-dlf fusion : 

2  o 

ne  >:>  electrons/cm-3  (15a) 


For  f ree-dif fusion : 


ng  <<:  48T/£^  electrons/cm^  (15b) 

From  the  above  criteria,  it  is  seen  that  for  a  gradient 
length  scale  of  the  order  of  1  cm  in  a  wake  plasma  up  to  a 
few  thousand  °K  temperature,  ambipolar  diffusion  will  pre- 
vail  at  electron  number  density  much  greater  .than  lOVcm  . 

5.  APPLICATION  TO  HYPERSONIC  FAR-WAKE 

The  near-wake  of  hypersonic  objects  is  generally 
dominated  by  a  high  Mach  number  (relative  to  the  ambient 
air),  anisotropic  compressible  flow  field  with  strong  temper¬ 
ature  and  density  gradients,  as  well  as  rapid  chemical  re- 
20  21 

actions.  ’  To  such  a  flow  field,  classical  turbulence 
theories  are  clearly  not  applicable. 

At  sufficiently  high  flight  Reynolds  number 
Reoo  5  Uoorr/voo  far  above  that  required  for  laminar- 

turbulent  transition  in  the  wake)  and  at  sufficiently  large 
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distances  downstream  of  the  object  (typically  a  few  hundred 
times  the  characteristic  dimension  r^  of  the  object),  how- 
ever,  ballistic  range  experiments  showed  that  the  gas 
density  fluctuation  as  inferred  from  the  film  contrast  of 

schlieren  photographs  tended  to  become  statistically  homo- 

* 

geneous  and  isotropic  within  the  turbulent  wake  volume. 

Thus,  the  possibility  exists  that  classical  isotropic  turbu¬ 
lence  theories  remain  applicable  to  at  least  some  part  of 
the  hypersonic  far-wake .  The  range  of  applicability  for  the 
classical  theories  will  then  be  limited  only  by  the  restric¬ 
tion  that  the  fluctuation  amplitudes  for  gas  temperature, 
density,  kinematic  viscosity,  and  scalar  diffusivity  all  be 
very  small  in  comparison  with  their  respective  mean  values, 
as  assumed  in  all  classical  turbulence  theories. 

In  turbulent  shear  flows,  the  fluctuation  ampli¬ 
tude  of  scalar  quantities  is  generally  of  the  same  order  of 
magnitude  as  the  averaged  "defect"  or  "excess"  of  such  quan¬ 
tities  within  the  turbulent  zone  from  the  ambient  condition. 


2l\ 

As  indicated  in  the  paper  of  Herrmann  et  al>  the  film  con¬ 
trast  correlation  actually  sampled  the  transverse  curvature 
of  the  refractive  index  distribution  along  the  averaged  ray 
path  within  the  turbulent  wake  volume,  and  hence  tended  to 
be  heavily  biased  toward  the  smaller  eddies  as  long  as  the 
density  fluctuation  associated  with  these  eddies  remained 
strong  enough  to  be  seen  by  the  schlieren  system.  Thus, 
the  observed  homogeneity  and  isotropy  in  density  fluctuation 
could  well  be  restricted  only  to  the  smaller  scale  structure 
of  the  wake  turbulence,  as  one  would  expect. 
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In  the  case  of  hypersonic  wakes,  the  temperature  excess 
AT  =  T  -  is  expected  to  remain  greater  than  or  comparable 
to  the  ambient  air  temperature  for  a  very  long  distance 
behind  the  object.  As  an  example,  in  the  case  of  quasi¬ 
equilibrium  flow  behind  a  hypersonic  sphere  of  22,000  ft/sec 
21 

velocity,  the  mean  wake  temperature  at  a  thousand  sphere 
radii  behind  the  sphere  (x/r  =  1*000)  is  expected  to  be 
approximately  1,500°K,  so  that  the  temperature  excess  would 
be  4  or  5  times  the  ambient  air  temperature  (depending  on 
altitude).  In  the  case  of  non-equilibrium  flow  (frozen 
dissociation),  the  temperature  excess  would  be  smaller,  but 
it  generally  would  still  take  a  distance  of  a  few  thousand 
sphere  radii  for  AT  to  subside  to  values  smaller  than  T  . 

Since  the  gas  density  p  is  roughly  inversely  pro¬ 
portional  to  the  gas  temperature  (for  constant-pressure  wake 

flow)  while  the  kinematic  viscosity  v  is  roughly  proportional 
2 

to  T  ,  the  condition  of  small  fluctuation  amplitudes  for  T, 
p,  and  v  does  not  appear  to  be  an  easy  condition  to  satisfy 
except  at  the  "very-f ar-wake" .  On  the  other  hand,  in  the 
very-far-wake  where  these  fluctuation  amplitudes  become  small, 
one  may  very  well  run  into  the  difficulty  of  transition  from 
ambipolar-dif fusion  to  free-dif fusion  for  the  electrons  as  ng 
decreases  [equations  (15a)  and  (15b)].  When  such  transition 
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occurs,  the  electron  diffusivity  D  becomes  a  rapidly  vary¬ 
ing  function  of  the  local  electron  density  ng(r). 

Prom  the  preceding  discussion,  it  is  quite  clear 
that  existing  turbulent  scalar  mixing  theories  may  find  very 
limited  range  of  applicable  conditions  in  hypersonic  wake 
plasmas  on  account  of  the  many  conflicting  restrictions. 
Nevertheless,  in  the  event  that  ambipolar-diffusion  does 
prevail  in  a  recombination-controlled  far-wake  (or,  alter¬ 
natively,  in  a  very-far-wake  where  the  dielectric  constant 
fluctuation  is  dominated  by  molecular  ion  density  fluctua¬ 
tions11),  then  one  shall  not  be  too  surprised  to  find  a 
universal  equilibrium  range  of  plasma  dielectric  constant 
fluctuation  with  a  spectral  shape  identical  to  that  observed 
in  low  speed  turbulence.  As  illustrated  in  Fig.  4,  the  shap 
of  such  a  spectrum  for  any  scalar  quantity  with  nearly  unity 
Prandtl  number  (i.e.,  v/D  =  1,  such  as  diffusion  of  heat  in 
perfect  gases  under  the  condition  of  small  temperature  fluc¬ 
tuation  amplitudes  or  ambipolar-diffusion  of  electrons)  is 

indeed  very  similar  to  the  turbulence  energy  spectrum,  as 
6  7  8 

Oboukhov,  Corrsin,  and  Batchelor0  all  predicted.  (See 
Section  III  above,  and  compare  Figs.  2  and  .) 

If  one  divides  the  scalar  fluctuation  spectrum 
shown  in  Fig.  4  for  Pr  =  0.7  by  the  square  of  the  wavenumber 
q  in  accordance  with  equation  (13),  one  then  obtains  the 
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universal  equilibrium  form  of  the  scattering  spectrum 

F(r,q).  If  the  spectral  separation  between  the  energy- 

containing  wavenumber  q^  and  the  Kolmogorov  wavenumber  q^ 

(which  is  roughly  the  same  as  the  Oboukhov-Corrsin  wave- 

\ 

number  q  in  this  case)  is  sufficiently  large,  such  scat- 

CJ  U 

tering  spectrum  would  be  characterized  by  a  power-law 
_  1  y 

segment  F  «  q“  '* corresponding  to  the  inertial-convection 

_S/ 

subrange  (r  «  q_/3)of  the  scalar  fluctuation  spectrum  up  to 

a  wavenumber  of  approximately  0.2qK,  and  dissipative  region 

beyond  q/q^  -  0.2  where  the  spectrum  falls  off  rapidly  with 

increasing  value  of  q/q^  at  an  ever  accelerating  pace.  As 

indicated  symbolically  in  Fig.  1,  the  existence  of  the 

inertial  subrange  in  the  energy  spectrum  is  conditioned  on 

a  sufficiently  high  value  of  the  turbulence  Reynolds  number 

(in  comparison  with  unity),  and  so  would  be  the  inertial- 

convection  subrange  in  the  scalar  fluctuation  spectrum  and 
_ii/ 

the  q  /J  segment  of  the  scattering  spectrum. 

The  preceding  discussion  on  the  spectral  form  of 
F(r,q)  in  the  universal  equilibrium  range  applies  only  to  a 
limited  spatial  region  about  the  point  £  within  the  wake. 
Since  the  radar  cross-section  a(0,XQ)  as  defined  in  equation 
(1)  is  given  by  an  integral  over  the  scattering  volume  V, 
the  wavelength-dependence  (or  frequency-dependence)  of  the 
radar  cross-section  cannot  be  used  directly  to  infer  the 
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spectral  form  of  F(r,q)  unless  the  spatial  variations  of 

— 5 - 

both  Ae  (r)  and  F(r,q)  happened  to  be  negligibly  small  ove 
the  scattering  volume. 

VI.  ESTIMATE  OF  THE  KOLMOGOROV  WAVENUMBER 

Within  the  universal  equilibrium  range  discussed 
in  the  preceding  section,  the  form  of  the  scattering  spec¬ 
trum  F(r,q)  depends  only  on  the  normalized  wavenumber  q/qv. 

It  also  follows  that  the  range  of  radar  wavelengths  over 
which  the  scattering  intensity  would  experience  the  effect 
of  diffusive  cutoff  to  the  scattering  spectrum  is  completely 
fixed  by  the  averaged  value  of  the  Kolmogorov  or  Oboukhov- 

O  l  O  V 

Corrsin  wavenumber  qK  =  (£/vJ)  -  (£/DJ)  within  the  scat¬ 

tering  volume  of  interest  at  any  given  scattering  angle  0 
(assuming  that  the  spatial  variation  of  F  is  sufficiently 
mild  throughout  the  scattering  volume). 

In  classical  turbulence  theory,  the  local  rate  of 
dissipation  of  turbulence  kinetic  energy  £  is  assumed  to  be 
a  given  quantity  characteristic  of  the  macroscopic  properties 
of  the  turbulence  field.  It  is  also  a  directly  measurable 
quantity  in  low-speed  turbulence  experiments.^  In  the  case 
of  hypersonic  wake  flow,  the  dissipation  rate  £  has  never 
been  directly  measured,  and  has  only  been  crudely  estimated 
on  occasion.  One  such  estimate,  for  the  case  cf  hypersonic 
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spheres,  was  reported  in  Ref.  26.  This  estimate,  which  was 
based  on  a  simple  extrapolation  of  low-speed  turbulence  data 
and  on  experimentally  observed  wake  growth  law  of  hypersonic 
spheres  in  ballistic  ranges,  gave  for  the  dissipation  rate 
in  the  far  wake 


£  *  0.06(x/rn)"2U^/rn  (16) 

where  x  is  the  distance  behind  the  sphere  along  the  wake 
axis,  rn  is  the  radius  of  the  sphere,  and  is  the  velocity 
of  the  sphere  relative  to  the  ambient  air.  It  is  interesting 
to  note  that  the  dissipation  rate  so  estimated  turned  out  to 
be  not  very  different  from  that  measured  by  Gibson  et  al1^ 
in  a  low  speed  wake  behind  a  sphere  in  a  water  tunnel,  which 
gave,  for  the  limited  range  of  value  of  x/r^  covered  by  this 
latter  experiment, 

£  =  0.27(x/rn)"2'35  U3//rn  .  (17) 

By  substituting  (16)  into  equation  (3),  and  defining  a 
"flight  Reynolds  number"  for  the  sphere  Re^  =  U^r  /v^,  one 
obtains  the  following  estimate  for  the  Kolmogorov  wavenumber 
in  the  wake  of  a  hypersonic  sphere. 
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q„  *  0.5{x/r)~h  (vyv)V'1  Re  ^  (18) 

In  the  hypersonic  wake,  the  kinematic  viscosity  v  varies 

roughly  with  the  square  of  the  gas  temperature  T,  so  that  the 

above  estimate  of  the  Kolmogorov  wavenumber  is  meaningful 

only  if  the  temperature  fluctuation  amplitude  is  indeed  small 

in  comparison  with  the  local  mean  temperature  T  as  discussed 

}  —  \ 

earlier,  and  in  such  case,  one  may  let  (vB/v)  4  =  (T^/T)  2. 

The  mean  gas  temperature  in  a  reacting  turbulent  flow  field 
is  generally  a  very  difficult  quantity  to  calculate.  It 
depends  on  the  turbulent  transport  and  mixing  rates  as  well 
as  on  the  chemical  reaction  model  for  the  turbulent  flow 
field.  At  any  rate,  such  mean  temperature,  and  subsequently 
the  Kolmogorov  wavenumber,  have  been  calculated  according 
to  the  reacting  turbulent  flow  model  proposed  in  Ref.  26  for 
some  typical  hypersonic  flight  conditions.  The  result  is 
reproduced  here  in  Fig.  5.  It  is  interesting  to  note  that 
the  Kolmogorov  wavenumber  so  calculated  (dotted  curves)  turned 
out  to  be  a  somewhat  wavy  function  of  distance  behind  the  sphere 
instead  of  a  monotonically  decreasing  function.  This  was  due 
to  the  competing  effect  between  the  rate  of  decrease  of  the 
estimated  turbulence  energy  dissipation  rate,  and  the  rate  of 
decrease  of  kinematic  viscosity  accompanying  the  continuously 
decreasing  mean  wake  temperature. 
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Ignoring  the  detailed  variation  of  the  dotted 
curves  with  distance,  the  results  illustrated  in  Pig.  5 
indicate  that  the  Kolmogorov  wavenumber  in  the  far  wake  of 

_  o  _  I  V 

the  sphere  is  of  the  order  of  5*10  r  Re  \  Thus,  for  a 

n  00 

sphere  of  1  ft  radius  (r  =  30.5  cm)  at  22,000  ft/sec 
velocity  (^7  km/sec)  and  150,000  ft  altitude  (^50  km),  one 
obtains  -  0.9  cm"^.  For  radar  back  scattering,  8  =  tt , 
and  q  =  (^Tr/Xo)sin(7r/2 )  =  47t/Xq.  Since  the  dissipative 
range  of  the  scalar  spectrum  begins  at  about  q/q^  =  0.2 
(see  Figs.  2  and  4),  one  may  conclude  that  diffusive  cutoff 
to  the  scattering  spectrum  will  be  felt  at  radar  wavelength 
XQ  shorter  than  about  50  cm  according  to  this  estimate. 
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FIG.  2  EXPERIMENTALLY  OBSERVED  TURBULENCE  ENERGY 
SPECTRUM  IN  UNIVERSAL  EQUILIBRIUM"  RANGE. 
(REPRODUCED  FROM  FIG.  10  OF  REF.  16) 
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FIG.  3  ENERGY  SPECTRUM  AT  DIFFERENT  STAGES  OF  DECAY  OBSERVED 
IN  GRID -TURBULENCE  EXPERIMENT  OF  STEWARD  AND  TOWNSEND. 
x/M  IS  THE  DISTANCE  DOWNSTREAM  OF  THE  GRID  IN  UNITS  OF 
THE  GRID  SPACING  M.( REPRODUCED  FROM  FIG.  7.6  OF  REF  12.) 
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FIG.  4  COMPARISON  OF  SCALAR  FLUCTUATION  SPECTRA  IN  UNIVERSAL 
EQUILIBRIUM  RANGE  OBSERVED  IN  WIND  TUNNEL  AND  IN  WATER 
TUNNEL  FOR  SCALARS  OF  DIFFERENT  PRANDTL  NUMBER  Pr  ■  v/D. 
(REPRODUCED  FROM  FIG.  2  OF  REF.  25). 


FIG.  5.  ESTIMATED  KOLMOGOROV  WAVENUMBER  IN  THE 
HYPERSONIC  WAKE  OF  A  SPHERE  OF  I  FT.  RADIUS 
AS  A  FUNCTION  OF  ALTITUDE  AND  NORMALIZED 
AXIAL  DISTANCE. (REPRODUCED  FROM  FIG.7  0F 
REF  26.) 


